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PREFACE 


This  report  is  the  result  of  a  Phase  I  SBIR  feasibility  study  into  optimum  navigation  for 
an  extended  range  intercept  missile.  The  problem  is  to  obtain  an  optimal  guidance  law 
for  both  the  midcourse  and  the  terminal  guidance  phase,  where  the  different  constraints 
for  the  two  phases  are  satisfied.  Direct  application  of  the  optimal  guidance  theory  will 
result  in  a  nonlinear  two-point  boundary  value  problem,  and  the  complications  due  to 
the  constraints  make  an  analytical  development  of  a  guidance  and  control  law  usually 
impossible. 

The  conventional  approach  to  this  problem  is  to  depend  on  iterative  methods  in  the 
implementation  of  the  optimum  guidance  law  that  satisfies  the  constraints.  In  this 
report,  a  computational  scheme  has  been  developed  that  does  not  require  iteration,  and 
that  is  generally  implementable  as  an  on-board,  real-time  guidance  law.  In  this  scheme 
the  nonlinear  equations  are  replaced  by  first  order  linear  approximations,  with  a  set  of 
corrections  added  to  account  for  the  nonlinearities.  These  nonlinearities  are  made  to  be 
known  functions  of  time  by  using  the  information  from  the  inertial  reference  unit  (IRU). 


Two  low-order  examples  have  been  included  in  the  report,  in  order  to  demonstrate 
the  computational  technique.  The  navigation  law  performance  and  the  sensitivity  to 
errors  and  updates  for  both  the  midcourse  and  the  terminal  guidance  modes  have  been 
evaluated  using  ATS’s  adjoint  simulation. 


1.0  INTRODUCTION 


The  guidance  problem  for  an  extended  range  intercept  missile  is  to  obtain  an  optimal 
guidance  law  for  both  midcourse  and  terminal  guidance  phases.  These  missiles  typically 
include  command  update,  inertial  midcourse  guidance  and  active  or  semiactive  terminal 
guidance. 

The  function  of  the  midcourse  guidance  law  is  to  minimize  energy  loss  and  bring  the 
heading  error  to  zero  at  handover.  Handover  is  the  change  from  midcourse  to  terminal 
guidance  and  occurs  following  target  acquisition  by  the  missile  seeker.  A  midcourse 
guidance  phase  by  itself  is  generally  not  fast  enough  to  consistently  achieve  a  miss 
distance  that  is  within  the  lethal  radius  of  the  warhead  and  a  terminal  guidance  mode  is 
necessary  to  achieve  the  required  miss-distance  results.  The  use  of  midcourse  guidance 
followed  by  a  relatively  short  period  of  terminal  homing  offers  a  significant  improvement 
in  firepower  and  missile  intercept  coverage  at  the  expense  of  the  addition  of  inertial 
sensors  and  their  initialization  and  data  links.  The  success  or  failure  of  the  terminal 
guidance  phase  is  affected  by  several  factors.  Seeker  acquisition  plays  an  important  role 
at  the  desired  time  of  handover.  The  heading  error  at  handover  is  an  essential  variable, 
affecting  the  terminal  miss-distance,  depending  on  the  missile-to-target  range,  missile 
speed,  autopilot  and  control  loop  dynamic  response,  guidance  filter  time  constants  and 
possible  radome  boresight  error  coupling.  As  the  missile  approaches  intercept,  speed 
and  maneuverability  can  become  critical  if  a  crossing  and/or  maneuvering  target  has 
to  be  intercepted,  since  the  interceptor  requires  from  three  to  five  times  the  maneuver 
capability  of  the  target.  Other  factors  affecting  terminal  guidance  accuracy  and  miss- 
distance  results  include  target  glint  and  severe  fades  in  the  target-reflected  signal  at 
critical  times  prior  to  intercept. 

A  combined  guidance  law  for  the  midcourse  and  terminal  guidance  phase  is  developed  in 
this  report  using  optimal  control  techniques.  However,  direct  application  of  the  optimal 
control  theory  will  result  in  a  nonlinear  two-point  boundary  value  problem.  The  problem 
can  be  further  complicated  by  lift,  thrust  and  control  constraints  forced  by  stuctural  and 
angle-of-attack  limits  and  the  direct  analytical  development  of  a  guidance  and  control 
law  is  usually  out  of  reach. 

Some  effective,  sophisticated,  numerical  optimization  methods  have  been  programmed, 
often  in  double  precision  arithmetic,  that  construct  a  single  optimum  guided  trajectory 
in  10  to  20  minutes.  These  techniques  will  not  be  applicable  for  many  years  to  real¬ 
time  control  of  missiles  with  compact,  light-weight  navigation  and  control  computers. 
Attempts  have  been  made  to  pre-compute  one  nominal  optimum  trajectory  and  then 
fit  thrust  angle  regimes  in-flight  to  the  pre-computed  optimal  trajectory.  This  kind  of 
implicit  solution  i«  useful  if  the  perturbations  in  the  initial  conditions  and  perturbations 
in  the  thrust  level  are  sufficiently  low.  But  if  the  mission  requires  many  different  powered 
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flight  maneuvers,  such  as  with  maneuvering  re-entry  vehicles,  the  solution  based  on  pre¬ 
computer  reference  trajectories  becomes  very  unattractive. 

In  this  report  two  more  easily  implementable  computational  schemes  have  been  devel¬ 
oped  that  allow  real-time  on-board  implementation.  The  first  uses  a  quasi-analytical 
approach  that  provides  a  solution  process  without  the  traditional  dependence  on  itera¬ 
tive  numerical  methods,  and  the  second  takes  advantage  of  the  knowledge  of  the  value 
of  the  nonlinear  term,  obtained  from  the  available  inertial  reference  unit. 

Since  the  intercept  missiles  have  to  counter  a  wide  variety  of  targets  under  extreme 
operating  conditions,  it  is  required  that  the  optimal  navigation  law  is  (1)  independent 
of  standard  conditions  and  (2)  nominal  trajectories  and  (3)  transitions  smoothly  from 
midcourse  to  terminal  navigation. 

The  first  two  requirements  can  be  met  by  seeking  the  solution  of  the  actual  differential 
equations  of  motion,  while  the  third  requirement  has  been  met  by  letting  the  solution 
be  subject  to  the  appropriate  boundary  conditions  at  terminal  handover.  A  requirement 
for  minimum  energy  during  flight  can  be  included  in  a  performance  index  that  defines 
the  parameters  or  states  to  be  minimized.  Since  the  minimization  of  the  desired  pa¬ 
rameters  can  often  lead  to  conflicting  strategies  (i.e.  minimum  fuel  and  minimum  time) 
the  weighting  terms  will  be  employed  in  the  performance  index  (or  cost  functional)  to 
define  the  degree  of  minimization  for  each  parameter.  Ideally,  the  determination  of  the 
parameters  to  be  minimized  is  a  direct  result  of  the  mission  requirements.  In  certain 
cases  however,  the  translation  of  mission  requirements  to  the  performance  index  leads 
to  unrealizable  or  non-unique  navigation  strategies.  In  that  case,  additional  terms  (sec¬ 
ondary  mission  requirements)  have  to  be  added  to  the  performance  index.  Care  must 
be  taken  in  selecting  and  weighing  the  additional  performance  criteria,  such  that  the 
primary  criteria  are  not  significantly  de-emphasized. 

System/state  constraints  can  be  handled  explicitly  or  implicitly  by  the  optimal  control 
formulation.  In  the  vehicle  control  problem,  these  constraints  may  be  factors,  such 
as  acceleration  limits,  seeker  angle  or  rate  limits,  terminal  constraints,  such  as  zero 
miss-distance,  terminal  aspect  angle,  or  other  state  constraints  that  may  limit  the  per¬ 
formance  of  the  system.  If  those  constraints  are  included  explicitly,  the  solutions  lead 
to  a  two-point  boundary  value  problem  that  can  only  be  solved  by  iterative  numerical 
techniques  which  are  very  sensitive  to  modeling  and  measurement  error,  making  on-line 
real-time  solutions  with  microprocessors  difficult  at  best.  For  these  reasons,  explicit 
constraints  are  handled  implicitly  through  the  performance  criterion  by  minimizing  the 
constrained  factors. 

Two  low  order  examples  have  been  included  in  this  report  to  demonstrate  the  com¬ 
putational  schemes  used  to  allow  real-time  implementation.  Of  those  two,  the  quasi- 
analytical  solution  is  the  more  general  and  powerful  approach  and  although  the  tech- 


nique  may  not  be  applicable  for  point-to-point  guidance  techniques,  it  warrants  further 
research  for  application  to  intercept  missile  guidance  techniques. 

The  navigation  law  performance  and  sensitivity  to  errors  and  updates  have  been  eval¬ 
uated  using  ATS’s  adjoint  simulation  for  both  the  midcourse  and  terminal  guidance 
modes. 


2.0  OPTIMAL  MIDCOURSE  NAVIGATION  LAW  DEVELOPMENT 

2.1  Introduction 

As  stated  in  Section  1.0,  there  are  many  considerations  in  developing  the  optimal  navi¬ 
gation  laws  via  optimal  control  theory.  The  criteria  to  be  considered  include  complexity, 
ease  of  implementation,  sensitivity  to  errors  and  updates,  availability  of  measurements, 
ease  of  handover  to  terminal  navigation  for  the  midcourse  navigation  law,  and  mini¬ 
mum  miss  distance  for  the  terminal  navigation  law.  Navigation  law  performance  and 
sensitivity  to  errors  and  updates  have  been  evaluated  for  the  BIM  using  ATS’s  adjoint 
simulation  for  both  the  midcourse  and  terminal  navigation  modes.  For  the  terminal 
navigation  mode  the  resulting  miss  is  also  a  function  of  statistical  disturbances,  such 
as  glint  or  angular  scintillation  noise,  (which  increases  with  decreasing  range),  effective 
receiver  noise,  (which  decreases  with  decreasing  range),  and  amplitude  or  range  inde¬ 
pendent  noise  caused  by  the  servo  system  and  by  amplitude  fluctuations  in  the  received 
signal.  In  addition,  for  a  system  such  as  the  BIM,  parasitic  coupling  with  the  vehicle’s 
airframe  caused  by  aberration  of  the  electromagnetic  energy  as  it  passes  through  the 
cover  protecting  the  terminal  sensor,  needs  to  be  considered.  This  error  typically  de¬ 
pends  on  the  vehicle’s  altitude  and  it  couples  the  airframe  dynamics  to  the  terminal 
sensor  measurements,  causing  possible  airframe  instability.  Because  of  the  interest  ex¬ 
pressed  in  this  error  source  during  a  visit  to  BMO,  a  brief  analysis  demonstrating  the 
resulting  feedback  paths  that  can  cause  missile  system  instability  has  been  included. 

Historically,  miss  distance  simulations  for  a  system  with  statistical  disturbances  have 
employed  tedious  Monte  Carlo  techniques  to  obtain  terminal  miss-distance  results.  ATS 
has  used  its  adjoint  simulation  to  determine  the  effects  of  initial  condition  errors  at  han¬ 
dover  and  the  effects  of  statistical  disturbances  on  missile  miss-distance  performance. 
While  Monte  Carlo  simulations  require  many  runs  to  obtain  statistically  meaningful  re¬ 
sults,  the  adjoint  technique  of  analysis  is  a  proven  statistical  analysis  tool  that  generates 
the  rms  miss-distance  in  only  one  run.  The  adjoint  technique  of  analysis  and  the  adjoint 
simulation  used  in  this  Phase  I  Study  are  described  in  Appendix  A  .and  B. 

2.2  Midcourse  Navigation  Law  Formulation 

The  development  of  the  optimal  midcourse  navigation  policy  has  been  accomplished 
using  the  differential  rquations  of  motion  of  the  vehicle  in  the  navigation  frame,  that 
is  used  by  the  inertial  platform  and  was  transferred  at  launch  to  the  inertial  reference 
system,  together  with  the  handover  to  terminal  navigation,  which  in  turn  is  a  function 
of  the  target  location.  If  the  target  is  a  moving,  non-cooperative  target,  this  handover 
to  terminal  navigation  point  may  have  to  be  updated  during  its  midcourse  flight. 

The  derivations  assume  that  the  navigation  takes  place  in  an  earth-centered  inertial 
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cartesian  coordinate  reference  frame. 


The  midcourse  navigation  problem  is  a  two-point  boundary  value  problem  which  consists 
of  the  vehicle’s  current  position  and  velocity  vectors 

r{to)  =  [X01  2/o>  Zo)  (l) 

v(t0)  =  (x0,y0,  z0) 

and  the  required  position  and  velocity  factors  at  handover 

r(tf)  =  (zf,yf,zf)  (2) 

v(tf)  =  (xf’Vf’Zf) 

The  requirement  is  then  to  find  the  thrust  acceleration  program  for  t0  <  t  <  tf  that 
drives  the  vehicle  from  its  initial  boundary  condition  to  the  terminal  boundary  condition. 
The  differential  equations  of  motion  of  a  missile  are 

r  =  g  +  a  (3) 


or  in  component  form 


x  =  gx  +  ax 

y  =  gv  +  av  (4) 

z  =  gz  +  ax 


If  the  gravity  field  is  assumed  to  be  spherical,  the  above  equations  are  non-linear,  i.e. 


[IX 

-I-,  rt 

JU  — 

(x2  +  y2  +  *2)5 

V  U' 

4  /  — — 

My 

- 1  ft 

y  — 

(x2  +  y1  +  22)’ 

I  Uj 

y  — 

f LZ 

-  (  ft 

6  — 

(x2  +  y2  4-  22)5 

•  Ur. 

(5) 


The  difficulty  in  obtaining  an  optimal  control  solution  to  the  above  dynamic  equations  is 
the  presence  of  the  non-linear  terms.  Only  rarely  is  it  feasible  to  solve  the  resulting  two- 
point  boundary  value  differential  equation  for  a  non-linear  system,  and  the  development 
of  “exact”  guidance  and  control  schemes  for  non-linear  systems  is  usually  out  of  reach. 
In  this  report,  however,  we  will  obtain  the  optimal  control  solution  using: 

1.  a  quasi-analytical  approach  obtained  from  a  straightforward  expansion  of  pertur¬ 
bation  methods,  which  provides  a  solution  process  without  the  traditional  depen¬ 
dence  on  iterative  numerical  methods,  and 
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2.  a  guidance  and  control  scheme  where  the  guidance  computer  generates  the  non¬ 
linear  terms  as  a  function  of  missile  position,  and  allows  the  optimal  control  prob¬ 
lem  to  be  treated  as  a  linear-quadratic  system. 

For  both  methods  the  performance  criterion  used  is  minimum  mean-squared  error.  There 
are  a  number  of  reasons  for  this  choice.  As  is  always  true  in  optimum  linear  theory 
-  whether  optimum  filtering,  optimum  control  or  other  special  case  -  the  quadratic 
error  criterion,  of  all  reasonable  criteria,  leads  to  the  most  tractable  formulation  of 
the  optimum  design  problem.  In  addition,  a  criterion  is  desired  which  is  universally 
applicable  to  all  forms  of  input  signals.  No  other  such  criterion  has  demonstrated  to 
give  superior  results  over  a  broad  range  of  problems. 

2.3  Quasi- Analytical  Solution 

The  non-linear  equations  from  (5)  can  be  written  in  the  state  variable  form  as  follows: 


*i 
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or  in  matrix  notation 
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(7) 
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(») 


We  seek  the  optimal  control  trajectory  that  will  minimize  the  quadratic  performance 
index 

J  =  ^ xTSx\t=t/  +  1}  Jt  * \xtQx  +  uTRu)dt  (10) 

where  R  and  S  are  positive  definite  weighting  matrices  and  Q  a  symmetric  positive  semi- 
definite  weighting  matrix  in  which  Q  =  0  is  not  excluded.  The  Hamiltonian  formed  from 
the  system  equation  (7)  and  the  integrand  of  equation  (10)  is  given  by 

H  =  ^( xtQx  +  u  tRu)  +  XT{Ax  +  Bu  +  f{x))  (ll) 

where  the  costates  A  are  a  set  of  as  yet  undetermined  Lagrange  multipliers.  Pontryagin’s 
necessary  conditions  for  determining  the  optimal  control  obtained  by  operating  on  the 
Hamiltonian,  yield  the  equations 

x  =  H\  =  Ax  +  Bu  +  f{x) 

A  =  -H,  =  -Qx-ATX-[^-)TX 

0  =  Hu  —  Ru  +  Bt\ 

with  the  boundary  condition  A (tf)  =  Sx(tj). 

Solving  equation  (14)  for  u  results  in  u  =  —R~1BTX  and  substituting  this  into  equation 
(12)  reduces  the  optimal  control  problem  to  two-coupled  first  order,  non-linear,  ordinary 
differential  equations,  i.e. 

x  —  Ax  -  BR~1BtX  +  f(x)  (15) 

A  =  -Qx  -  At A  -  (^l)rA  (16) 

Combining  the  unknowns  x  and  A  into  a  single  augmented  state/costate  vector  z,  the 
optimal  control  problem  can  be  restated  as 


(12) 

(13) 

(14) 


z  =  Fz  +  ep 


(17) 


where 


F  = 


z  =  (XT  xT)T 

A  -BR~lBT 
-Q  -AT 


(18) 

(19) 


._(  /(*)  'l  ,  , 

"  -  (  -(^)7A  j  (20) 

with  boundary  conditions  x(t0)  =  i(0)  and  A (tf)  =  Sx(tj )  and  where  the  dimensionless 
parameter  €  is  a  bookkeeping  term  to  keep  track  of  the  order  of  the  nonlinear  terms. 
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After  solving  the  two-point  boundary  value  problem  given  by  equation  (17),  the  state 
trajectories  and  optimal  control  can  be  determined.  However,  because  of  the  presence  of 
the  non-linear  terms,  the  system  governed  by  equation  (17)  is  analytically  intractable. 
Although  there  are  many  iterative  techniques  available  for  solving  such  non-linear  sys¬ 
tems,  they  are  not  suitable  for  on-board  implementation.  We  will  therefore  develop  a 
quasi-analytical  technique  that  will  circumvent  the  iterative  techniques. 

Assume  that  the  solution  to  equation  (17)  may  be  represented  by  a  power  series  in  terms 
of  a  small  parameter  e  by 

z{t)  =  z0(t )  -|-  ezi{t)  +  e2z2(t)  +  e3z3(f)  +  ...  (21) 

For  small  nonlinearities  the  series  will  converge,  while  the  accuracy  will  improve  as 
the  nonlinearities  approach  zero.  Similarly,  for  a  convergent  series, the  accuracy  can  be 
improved  by  using  more  terms,  in  fact  as  the  number  of  terms  in  the  series  approaches 
infinity,  the  solution  given  by  equation  (21)  will  be  exact.  Substituting  equation  (21) 
into  (17)  we  can  write 

z0  +  ezi  4-  e2z2  +  e3(0)  =  F  z0  +  eF  z2  +  e2Fz2  +  £Pi(z0)  +  e2p2(z0,z1 )  +  e3(0)  (22) 

where  the  nonlinear  terms  have  been  expanded  in  a  similar  power  series  and  the  depen¬ 
dence  of  the  nonlinear  term  on  the  z  variable  is  indicated  by  z,  .  Equating  terms  with 
equivalent  powers  of  e  results  in  the  following  set  of  differential  equations. 


Zo 

=  Fz0 

(23) 

Zl 

=  Fzi  +  P\{z0) 

(24) 

Z2 

=  Fz2  +  p2{z0,z1) 

(25) 

with  the  boundary  conditions 


x(*f)  \ 

K(tf )  ) 


and 

Zi^  =  (  A,(f0)  )  =  (  0  )  *  =  1’2’3’- 

and  where  the  final  conditions  of  the  states  and  costates  are  related  to  each  other 
through  the  boundary  condition  A,(£/)  =  Sz^tf). 

Thus,  we  obtained  a  strictly  linear  first  order  approximation  (equation  (22))  to  the 
original  nonlinear  problem  and  a  series  of  correction  terms  (equations  (23)  through 
(25))  to  account  for  the  effects  of  the  nonlinearities,  while  the  nonlinear  term  in  the 
correction  terms  is  independent  of  the  z  variable  of  that  particular  differential  equation. 
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Furthermore,  the  order  can  simply  be  extended  to  achieve  the  degree  of  accuracy  required 
for  a  specific  problem. 

The  solution  of  the  system  of  equations  given  by  equations  (23)  through  (25)  is  given 
by  t 

Zi(t)  =  eF(t)z,(0)  +  f  ef(1-r)p,(r)dr  t  =  0,1,2,...  (28) 

J  o 

or  equivalently 

Zi{t)  =  $(t)*i( 0)  +  $(t  -  T)pi(r)dt  i  =  0,1,2,...  (29) 

with  p0(t)  =  0  and  where  without  loss  of  generality  t0  has  been  replaced  by  0.  In 
shorthand  notation  we  can  rewrite  (28)  as 

z{(t)  =  $(£)*<  (0)  +  Vi(t)pi  (30) 

However,  at  this  point  the  initial  costates  A,(f0)  thus  z<(£0)  are  as  yet  undetermined.  If 
we  expand  the  state  transition  matrix  into 


and  the  last  term  of  equation  (30)  into 


= ( *;:l) ) 


then  on  recalling  the  boundary  conditions  of  equation  (20)  we  can  write  at  t  —  tf 

xi(tf)  =  $ii(£/)®.(0)  +  $i2(£/)A,(0)  4-  Vu(tf)  (33) 

SXi(tf)  =  $2l(£/)a:t(0)  +  ^22  (£/)  A,'  (0)  +  ^2.(£/) 

in  which  A,(0)  is  the  only  unknown.  Multiplying  equation  (33)  by  the  positive  definite 
matrix  S  and  subtracting  the  resulting  equation  from  (34)  results  in 

[*22(*/)  -  S*„{tf)) A,(0)  =  (5»11(*/)  -  $2i(£/)]x,(0)  +  S9u{tf)  -  *»(*,)  (34) 

from  which  the  original  costates  are: 

\(0)  =:  [^22 (^/)  —  ,S,^i2 (^/)]  1[(*S'^ii(£/)  _  ^2i(£/))xt(0)  +  —  ^2, -(£/)]  (35) 

Now  with  all  initial  conditions  known,  equation  (29)  can  be  used  to  obtain  the  optimal 
control  at  any  time  in  the  interval  0  <£<£/.  The  original  nonlinear  control  problem 
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has  thus  been  expanded  to  a  set  of  non-homogeneous,  linear,  optimal  control  problems 
that  may  be  solved  sequentially. 

The  effectiveness  of  the  quasi-analytical  method  is  demonstrated  with  an  example  of  a 
system  with  both  quadratic  and  cubic  nonlinear  terms.  The  system  is  given  by 

ii  =  x2  (36) 

x2  =  —Xi  —  .lx2  +  u  +  .l(uXi  —  .5ZJ3)  (37) 

or 

-1.!  )  z  +  (  .1  )  tt  +  (  .i  )  (uzi  "  -5zi3)  (38) 

or  in  short 

x  =  Ax  +  bu  +  cf(x,  u)  (39) 

The  objective  is  to  determine  the  optimal  controls  that  will  drive  zx  from  one  to  zero  in 
a  2  second  interval.  The  corresponding  performance  criterion  is 

J  =  ^ xTSx+  ^  J  u2dt  (40) 

with  sx  =  s2  =  100.  The  Hamiltonian  is  given  by  H  =  \u2  4-  A T(Ax  +  bu  +  cf(x,  u))  and 
Pontryagin’s  necessary  conditions  for  determining  the  optimal  controls  are 


Hx 

=  x  =  Ax  +  bu  +  cf(x,u) 

(41) 

-Hx 

=  A  =  -At\  — 

(42) 

Hu 

—  0  =  u  +  bTX  +  w^-cTA 
du 

(43) 

or 

u  =  —bTX  —  -^-cT  A 
ou 

■ 

(44) 

Substituting 

the  optimal  control 

into  the  differential  equations  results  in: 

x  =  Ax  —  bbT  A  —  +  c/(z,  u) 

(45) 

A  =  - 

.AtX-c^ctx 

ox 

Combining  the  unknowns  x  and  A  into  a  single  augmented  state/costate  vector  z,  the 
optimal  control  problem  can  be  restated  as 

z  =  Fz  +  ep  (46) 


where 


z  =  (\T  xT)T 
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(48) 


0  10  0' 

-1  -.1  0  -1 

0  0  0  1 

0  0-1  -.1  J 

and 

0 

—  3^2  +  ( UXi  —  .5S!3) 

0 

-(u  —  1.5xi2)A2 

The  effectiveness  of  the  optimal  control  approximation  was  evaluated  by  integrating 
equation  (46)  numerically  using  a  fourth  order  Runga  Kutta  routine  and  examining  the 
terminal  boundary  condition  errors  of  the  numerically  integrated  solution  with  those  of 
a  second  order  expansion. 

The  second  order  expansion  yields  a  final  condition  error  of  x(t)  =  —.000321  from  the 
integrated  equation  of  motion.  Although  not  exactly  zero,  the  error  is  less  than  0.035%. 
By  comparison,  the  linearized  optimal  control,  obtained  by  dropping  the  nonlinear  terms 
produces  x(tf)  =  —.0382,  or  an  error  of  slightly  less  than  4.0%.  Thus  the  perturbation 
method  reduced  the  error  by  more  than  two  orders  of  magnitude  for  a  second  order 
expansion,  demonstrating  the  effectiveness  of  the  method. 

The  secret  to  the  success  of  this  powerful  method  lies  of  course  in  the  selection  of  a 
convergent  series  for  the  assumed  solution.  This  will  occur  automatically  if  the  selection 
of  z0  is  close  to  the  actual  solution,  ensuring  that  the  correction  terms  will  be  small. 
Selection  of  such  an  initial  solution  is  illustrated  by  an  example  of  determining  the 
near-circular  orbit  period  of  an  earth  satellite.  The  coupled  radial  and  tangential  force 
equations  for  such  a  system  are  given  by 

f  —  rip2  =  (50) 

and 

rxp  +  2rxjj  =  0  (51) 

Since  this  latter  equation  can  be  written  as  ^(r2^)  =  0,  we  can  write  r2^  =  h  where 
we  defined: 

r  as  the  distance  from  the  center  of  the  earth  to  the  satellite, 

as  the  angle  between  r  and  an  arbitrary  reference, 

H  as  a  constant  defining  the  specific  gravitational  force, 

h  as  the  constant  specific  angular  momentum  of  the  orbit. 
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Eliminating  xp  from  the  above  equations  results  in  the  non-linear  system  differential 
equation 

h 2  n 

T  r3  ~  r2 

for  which  we  have  to  assume  a  solution.  The  geometry  of  the  problem  suggests  a  solution 
of  the  form 

r  =  R  +  6r  cos  xp  (53) 

where  xp  =  uit  with  R  and  SR  as  constants. 

Substituting  this  expression  in  the  nonlinear  system  equation  results  in 

-  u26r  cos  *P  ~  ^3  (1  +  ^cosip)-3  =  - -£2  (i  +  —cosrp)-2  (54) 

The  assumption  of  a  near-circular  orbit  implies  that  SRj R  <§;  1  and  the  above  equation 
can  be  approximated  by  the  following  first-order  expansions 


In  a  missile  system  that  is  aware  of  its  physical  position  with  its  Inertial  Reference  Unit 
(IRU),  we  have  the  option  of  calculating  the  actual  value  of  the  nonlinear  terms  such  as 
gravity  at  each  instant  of  time,  and  use  that  value  to  modify  the  required  total  optimal 
acceleration.  This  total  optimal  acceleration  can  be  calculated  using  the  linear-quadratic 
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approach  as  discussed  below.  Using  the  total  acceleration  concept  we  can  rewrite  the 
system  equations  (5)  as: 


or  in  matrix  notation 


where 


(ij3  +  x32  -f  x62)3/2 

x4 

_ -/^3 

(Xi2  +  X32  +  Z52)3/2 
*6 

_ -M»5 

(xx2  +  X32  +  X52)3/2 


x  =  Ax  +  Bu 


+  az  =  Uj 


+  av  =  u2 


+  a,  =  u3 


0  1  0  0  0  0 
0  0  0  0  0  0 
0  0  0  1  0  0 
0  0  0  0  0  0 
0  0  0  0  0  1 
0  0  0  0  0  0 


and  B  = 


0  0  0 
1  0  0 
0  0  0 
0  10 
0  0  0 
0  0  1 


We  seek  the  optimal  control  that  will  minimize  the  quadratic  performance  index 


J  =  1/2 xT(tf)Sx(tj)  +  1/2^  * (uT Ru)dt 


where 


Si  0 

0  s2 


0  0 
s2  0 

0  83 


0  0 
0  0 


s  5  0 

0  5ft 


and  R  = 


are  positive  definite,  diagonal  weighting  matrices. 

The  Hamiltonian  formed  from  the  system  equations  (61)  and  the  integrand  of  the  equa¬ 


tion  (63)  is  given  by 


H  =  1/2 utRu  4-  \T [Ax  +  Bu) 


where  A  is  the  costate  vector  of  the  Lagrange  multipliers.  Pontryagin’s  necessary  con¬ 
ditions  for  determining  the  optimal  control  obtained  by  operating  on  the  Hamiltonian 
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yield  the  equations 


x  =  Hx  =  Ax  +  Bu  (66) 

A  =  -Hz  =  -At  A  (67) 

0  =  Hu  =  Ru  +  BT  A  (68) 

with  the  boundary  conditions  x(t0)  =  z(0)  and  A(t/)  =  Sx(tf).  Solving  equation  (68) 
for  u  results  in  u  =  —R~1BT A  and  substitution  in  equation  (66)  reduces  the  optimal 
control  problem  to  two-coupled  first  order,  ordinary  differential  equations,  i.e. 

x  =  Ax-BR~lBT  A  (69) 

A  =  -At  A  (70) 

with  boundary  conditions  x(t0 )  =  x(0)  and  A(t^)  =  5x(f/).  Writing  out  the  A  part  we 
have 

Ax  =  0 
A2  =  — Aj 

A3  =  0 

^4  =  -^3  (71) 

Ag  =  0 

A6  =  — A6 


with  boundary  conditions  A (tf)  —  Sjx(tf).  Integration  from  ts  to  t  gives 


Ai(0  =  «l*l(f/) 

A2(t)  =  51x1(t/)(t/  -t)  +s2x2(tf) 
A3(t)  =  33X3(tf) 

A 4(0  =  «3*3(</)(*/  -  t)  +  s4xA{tf) 
A$(t)  =  s6x3(tf) 

A e(t)  =  *b*5(*/)(*/  ~  0  +  s6x6(tf) 
and  the  optimal  control  using  u°  =  —R~1BT A  is  then 


«1°  =  -£xi(*/)(*/  -0  -  y?*a(*/) 

»2°  =  ~^*3(«/)(</  -0  “  £*<(*/) 

U30  =  “£*#(*/)(*/  -  0  -  £*e(*/) 


(72) 


(73) 


or 


u°(t)  = 


/  -ff  0  0 

0  0  -%T 

V  0  0  0  0 
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0 

0 


—  liT 

r» 


0 

0 


_£A 


\ 

*(*/) 

/ 


(74) 


which  gives  us  the  optimal  control  in  terms  of  the  state  variables  at  final  time  tj  and 
time-to-go,  T  —  (tf  —  t). 

To  obtain  the  optimal  control  in  terms  of  the  current  time  t,  we  use  the  differential 
equations  of  motion  which  can  be  written  as 

Xj  =  X3 

x2  =  ~Xl(tf){t/-t)’-^Xi(tf) 

x3  =  XA 

Xi  =  -^*3(*/)(</-0"*4(*/)  (75) 

is  =  x6 

*«  =  -^®b(*/)(*/  ~  0  -  ^x6(tf) 

with  boundary  conditions  at  t0  =  x(f0).  Integrating  the  above  equation  from  t0  to  tj 
results  in 

*1  (</)  =  +  +  *,(») 

*>(«/)  =  +  r  +  x,(o)  (76) 

*<(«/)  =  -^*3(f/)rI-|*<(i,)r  + *,(<>) 

*.(»/)  =  -^*s(‘/)r3-^x,((/)T1  +  x6(o)r  +  xt(o) 

*«(</)  =  + 

where  T  =  tf  —  ty  or  in  matrix  notation 

1  +  fcT3  ftT2  0  0  0  0 

ftT3  1  +  ffT  0  0  0  0 

0  0  1  +  ^r3  ^r2  0  o1 

0  0  ^T2  1  +  0  0 

0  0  0  01  +  ^r3  ^r2 

0000  ^r2  1  +  ^r2 

f  1  T  0  0  0  0  ' 

0  1  0  0  0  0 

0  0  1  T  0  0  (78) 

0  0  0  1  0  0  *' 

0  0  0  0  1  T 

0  0  0  0  0  1  , 
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x(tf)  (77) 


If  the  above  equation  can  be  solved  for  x(tf ),  then  the  result  can  be  substituted  into 
Eqn  (74)  to  obtain  the  optimal  control. 

By  rewriting  the  coefficient  matrix  of  Eqn  (78)  as 

'll  12  0  0  0  0  ^ 

21  22  0  0  0  0 

0  0  33  34  0  0 

0  0  43  44  0  0 

0  0  0  0  55  56 

^  0  0  0  0  65  66  j 

its  inverse  has  been  obtained  and  is 


( 

22 

X 

12 

~x 

0 

0 

0 

0 

\ 

21 

“X 

11 

X 

0 

0 

0 

0 

0 

0 

44 

34 

0 

0 

0 

0 

43 

-~E 

33 

0 

0 

0 

0 

0 

0 

66 

u 

56 

~77 

\ 

0 

0 

0 

0 

65 

55 

V 

where 

A  =  11x22  -  12x21 
B  =  33x44  -  34x43 
C  =  55x66  -  65x56 

Post-multiplying  Eqn  (74)  with  Eqn  (80)  and  the  right  hand  side  of  Eqn  (77)  or  Eqn 
(78),  substituting  the  actual  terms  from  the  coefficient  matrix,  and  letting  the  current 
time  t  =  t0,  results  in  the  continuous  optimal  feedback  control  terms 


Ul°(t)  = 


u2°(t)  = 


u3°(t)  = 


6  3lT(s2T  +  2rj) 

3i32T *  4  4r1s1T3  +  12rjS2r  4  12r^Xl[  > 
sl32T3  4  Ss^T2  -I-  Zrls2  _  /iX 
3232T*  +  4r1s1  T3  +  12  rxs2T  +  V2rx*X*W 
6 s3T{s4T  4  2r2)  _ 

s3s4T*  4  4 r2s3T3  4  12r2s,r  4  I2rfx*l) 

,  S3S4T3  4  3s3r2r2  -I-  3r2s4  _ 

s334T*  4  4r2s3T'3  4  12r2s,r  4  12 r2*Xi[t) 

_ 6 sbT{seT  4  2r3) _ 

35s6T*  4  4r33$T3  4  12 r3s6T  4  12 r32Zsl  j 
ssSqT3  4  3s5r3T 2  4-  3r3s6 
s5s6T*  4  4r3s5r3  4  12r3s6r  4  12r32Xeltj 
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(81) 


The  required  acceleration  along  each  axis  can  now  be  obtained  by  making  the  sum  of 
the  actual  and  gravity  acceleration  equal  to  the  total  required  acceleration,  thus 


u  °{t)  =  a(t)  +  g{t)  (82) 

or 

a(t)  =  u°(t)  -  g(t )  (83) 

for  each  component.  By  using  u^t)  —  gi(t);  u2(t)  —  gr2 (^ )  and  u3(t)  —  g$(t)  it  is  possible  to 
determine  the  optimal  acceleration  for  each  axis.  Furthermore,  constraining  equations 
such  as  the  constraining  relationship  between  the  length  of  a  and  the  three  components 
of  a  can  readily  be  added.  For  example,  if  an  engine  is  throttleable  as  well  as  gimballed, 
such  a  constraint  can  be  satisfied  by  ensuring  that 

Igt|  =  iaTx2  +  aTv2  +  aTz2) 5  (84) 

where  the  direction  cosines  of  the  desired  thrust  direction  given  by 

cos  O!  =  ~~~  (85) 

cos/?  =  (86) 

cosi=^f  (87) 

The  above  commands  have  been  calculated  in  the  inertial  frame  and  need  to  be  trans¬ 

formed  to  the  body  frame,  before  they  can  be  applied  to  the  air-frame,  as  discussed  in 
Section  3.5.  For  rt  =  r2  =  r3  =  r,  =  s3  =  =  s,  and  s2  =  s4  =  s6  =  s  the  general 

optimal  navigation  law  can  be  written  as 

oUX  6sT(sT  +  2r) 

U<  W  "  s2T*  +  4 rsT3  +  UrsT  +  12r»X<^  (88) 

4(a2T3  +  3  sT2  +  3  rs)  _ 
s2T*  +  4  rsT3  4-  12  rsT  +  12r>X<+1^* 

To  demonstrate  the  linear-quadratic  approach,  an  example  is  presented  of  a  low  order 
nonlinear  system  given  by 


x  =  v 

1 

v  =  a - ?  =  u 

x2 

with  the  performance  criterion 


J  = 


^xTSx  +  ^  J  u2dt 
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(89) 


(90) 


From  the  previous  section,  the  optimal  control  for  r  =  1  is  given  by 


“"(*)  = 


and  for  =  s2  =  2 


_ 6siT(s2T  +  2) _  ,  * 

slS2T*  +  4 5iT3  +  12s2T  +  12  U  J 
4(s1s2T3  +  3s1T2  +  3s2)  ,  , 

s1s2TH4s1Pf  12s2T+  12Z21  J 


u°(0  = 


6T(T  +  1) 

T*  +  2T3  +  6T  4-  3 
4T3  +  6T2  +  6 
T 4  +  2T3  -f  6T  +  3 


®i(0 

*2(0 


(91) 


(92) 


This  optimal  control  u  now  has  to  satisfy  the  condition  that  u°  =  —  ^  +  a,  i.e.  of  the 
total  control  requirement  u0,  part  —  is  provided  by  the  dynamics  of  the  vehicle,  while 
the  remaining  part  a  =  u  +  —■—  needs  to  be  provided  by  the  guidance  computer. 


2.5  Midcourse  Navigation  Law  Evaluation 

The  above  optimal  navigation  law  has  been  implemented  on  ATS’s  adjoint  simulation 
(See  Appendix  A  for  a  description  of  the  simulation).  A  second  order  autopilot  with 
damping  ration  of  .7  and  a  natural  frequency  of  5  rad/sec  and  an  airframe  time  constant 
of  2.5  sec,  was  used  in  the  studies. 

Miss  distance  sensitivity  runs  with  changes  in  the  location  of  the  handover-to-terminal 
navigation  point  caused  by  updates  in  the  point  from  the  ground  are  shown  in  Figure  1. 


The  results  show  that  with  no  update  or  with  an  update  prior  to  5  seconds-to-go  the 
resulting  miss  distance  is  negligible.  Worst  case  updates  would  occur  if  they  were  trans¬ 
mitted  at  .25,  .5  or  1.5  seconds-to-go. 

The  midcourse  navigation  law  development  and  the  above  miss-distance  results  assume 
that  there  are  no  errors  in  the  inertial  navigation  system  (INS)  of  the  vehicle.  In  real  life 
however,  errors  caused  by  initialization  or  instrument  drifts  of  the  INS  will  contribute 
to  the  miss  distance  at  handover.  Determination  of  those  errors  will  require  an  error 
analysis  of  the  INS,  analysis  of  the  transfer  alignment  from  the  master  inertial  reference 
system  to  the  INS  before  launch,  which  may  include  the  estimation  and  updating  of  the 
INS  gyro  drifts  and  accelerometer  biases.  However,  should  the  results  of  such  a  study 
and  the  mission  requirements  indicate  that  a  more  accurate  INS  is  required,  it  could 
be  achieved  by  updating  the  vehicle’s  position  from  the  ground  and  so  correct  for  any 
INS  errors  that  may  have  accrued  during  midcourse  flight.  The  timing  requirements  for 
this  INS  update  will  follow  a  similar  pattern  as  shown  in  Figure  1.  Figure  2  shows  the 
sensitivity  of  guidance  error  to  update  accuracy  as  a  function  of  range. 
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Fig.  1.  Miss-Distance  Sensitivity  with  a  Step  Displacement 
in  Handover-to-Terminal  Point  versus  Time-to-Go. 
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Fig.  2:  Guidance  Error  Sensitivity  to  Update  Accuracy 
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The  midcourse  guidance  errors  will  result  in  a  missile-to-target  line-of-sight  pointing 
error  during  the  acquisition  phase  at  handover  and  are  reflected  in  the  seeker  field  of 
view  requirements. 


If  R  is  the  missile-to-target  range  vector  at  launch  and  M  and  T  the  missile  and  target 
position  vectors  at  handover  to  terminal  guidance,  then  RMT  =  R  —  M  —  T  is  the  range 
vector  from  missile  to  target  at  handover.  If  at  handover  to  terminal  guidance  the  error 
in  target  position  due  to  the  target  prediction  error  is  AT  =  TA6t  and  if  the  error  in 
missile  location  at  that  time  is  AM  =  MA9m  then  the  maximum  look  angle  error  is 


_  AM+AT 
~  RMT 


MAOfrf  -t-  T  AO’f 

vJTo 


(93) 


and  if  A 9M  and  A 0T  are  independent  normal  random  variables,  then 


_2  _ 
°  A/9  — 


M\ 7  V  +  T2o 2 


OT 


V'Tgo2 

For  a  balanced  system  where  Mo6M  =  TagT  =  Rjoe  the  rms  look-angle  error  is 

\f2RjOg 

-  v,r8o 


(94) 


(95) 


where 


Ri  is  the  range  to  intercept 

Ve  is  the  missile-target  closing  velocity 

Tgo  is  the  time-to-go 


The  probability  that  the  radius  of  the  field-of-view  requirement  r  is  less  than  a  given 
value  R  is  defined  by 

P(r  <  iE)  =  1  -  e**/2***'  •  (96) 

or 

R  =  -y/2c&p[-ln(l  -  P(r  <  i?))]5  (97) 

Substituting  a ^  from  Eqn  (95)  gives 

OD 

R  =  v^Toc^-,nl-1-pl'r<R^i  (98> 

and  the  total  field  of  view  requirement  is  then 

D  =  'Sro°‘“{-ln{1-p(r<R))]i 
21 


(99) 


For  an  Rr  of  200,000  ft,  a  Vc  of  10,000  ft /sec,  and  a  oA6  of  10  millirad  the  field  of  view 
requirements  at  various  times-to-go  and  for  probabilities  of  target  containment  of  .9,  .99 
and  .999  are  shown  in  Table  I. 


Probability 
of  Containment 

Tgo 

sec 

FOV 

degrees 

.9 

20 

3.5 

10 

7.0 

5 

13.9 

.99 

20 

4.9 

10 

9.8 

5 

19.7 

.999 

20 

6.0 

10 

12.5 

5 

24.1 

Table  I:  Seeker  Field  of  View  Requirements 
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3.0  OPTIMAL  TERMINAL  NAVIGATION  LAW  DEVELOPMENT 


3.1  Introduction 

Development  of  the  terminal  navigation  law  closely  parallels  the  development  of  the 
midcourse  navigation  law.  However,  because  of  a  possible  moving  target,  the  differential 
equations  of  the  target  have  been  included  in  the  dynamic  model  equations,  increasing 
the  number  of  state  variables  from  six  to  nine.  In  addition,  since  the  final  time  is  now 
a  free  parameter,  it  has  to  be  estimated  and  methods  to  estimate  time-to-go  have  to 
be  developed.  Furthermore,  since  terminal  navigation  information  is  obtained  from  a 
terminal  sensor,  the  stochastic  aspects  of  this  sensor,  as  well  as  boresight  error  slope 
(discussed  below)  have  to  be  included  in  the  evaluation  of  the  resulting  navigation  law. 

3.2  Terminal  Navigation  Law  Formulation 

Similar  to  the  midcourse  navigation  law,  the  two-point  boundary  value  problem  for  the 
terminal  navigation  law  consists  of  the  vehicle’s  current  position  and  velocity  vectors 


r{t0)  =  {x0,y0,z0)  (100) 

v{t0)  =  (x0,y0,z0)  (101) 

and  the  required  position  and  velocity  vectors  at  intercept 

r(*/)  =  (*/,»/,*/)  (102) 

v{tf)  =  (xf,yf,Zf)  (103) 


The  requirement  is  to  find  the  acceleration  vector  program  for  t0  <  t  <tf  that  drives  the 
vehicle  from  its  present  position  to  intercept.  As  before,  the  total  required  acceleration 
profile  will  be  obtained  and  the  nonlinear  gravity  terms  are  then  taken  into  account  in 
a  way  that  necessitates  no  approximations. 

If  we  define  the  state  vector  with 

Zj  =  the  target /vehicle  relative  position  in  the  x  direction 

x2  =  the  target /vehicle  relative  position  in  the  y  direction 

z3  =  the  target /vehicle  relative  position  in  the  z  direction 

x4  =  the  target/vehicle  relative  velocity  in  the  x  direction 

z5  =  the  target/vehicle  relative  velocity  in  the  y  direction 

x6  =  the  target/vehicle  relative  velocity  in  the  z  direction 

x7  =  the  target/vehicle  relative  acceleration  in  the  x  direction 

x8  =  the  target/vehicle  relative  acceleration  in  the  y  direction 

x9  =  the  target /vehicle  relative  acceleration  in  the  z  direction 

and  if  we  further  assume  that  the  target  acceleration  can  be  modeled  as  a  first  order 
process,  i.e.  dj  =  —rax  then  we  can  formulate  a  set  of  first  order  linear  differential 
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equations  which  describe  the  system  dynamics  in  state  variable  form,  i.e. 


x\ 

= 

Xi 

x2 

= 

Xs 

Xs 

= 

Xs 

X< 

= 

x7  —  a. 

Xs 

= 

x8  —  a, 

x$ 

= 

Xq  “  flj 

x 7 

= 

—  TiX7 

Xs 

= 

—  T2Xg 

X> 

= 

T3®9 

(104) 


If  the  control  vector,  u,  is  defined  to  be  the  missile  acceleration,  then  the  state  model 
can  be  written  such  that 


where 


and 


x  =  Ax  +  Bu 


A  = 


O3  J3  O3  ^ 

(  °3  \ 

O3  O3  I3 

B  = 

~h 

O3  O3  -T\h  j 

l  03  J 

x 


f  xi  > 

x2 

*3 

*4 

xi 

x6 

x7 

x8 

Xo 


u 


az 

av 

a. 


(105) 

(106) 


(107) 


3.3  Optimal  Control  Formulation 

The  performance  criterion  to  be  minized  is  given  by 

J  =  ^ xT{tf)Sfx(tf )  +  \Jt  * (uTRu)dt 

where 


( 

«3 

0 

0  ^ 

( 31 

0 

0  > 

f 

ri 

0 

0  y 

s  = 

0 

0 

0 

with  S3  =  I  0 

s2 

0 

and 

R  = 

0 

0 

0 

0 

0  ) 

l  0 

0 

*3  ) 

0 

0 

r3 ; 
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(108) 


(109) 


This  cost  functional  has  no  weighting  on  the  final  relative  velocity  nor  on  the  target 
acceleration,  but  has  a  weighted  cost  on  the  control  or  vehicle  acceleration  through  the 
integral  term.  In  this  Phase  I  Study,  each  of  the  weighting  parameters  will  be  treated 
as  a  generic  variable  and  will  thus  appear  explicitly  in  the  solution. 

3.4  Optimal  Control  Solution 

Using  Eqns  (105)  and  (108)  the  Hamiltonian  is  constructed  as  follows 

H  =  ^uT  Ru  +  XT(Ax  +  Bu)  (HO) 

where  A  is  now  a  co-state  vector  with  dimension  9x1.  The  necessary  conditions  for 
optimality  are  then 

A  =  -Hx  =  -At  A  (111) 

0  =  Hu  =  Ru  +  Bt  A  (112) 

This  last  equation  can  be  re-written  as 

u°  =  —R~1BtX  (113) 

Substituting  Eqn  (113)  into  Eqn  (105)  provides 

x  =  Ax  -  BR~lBTX  (114) 

Thus  we  can  write  Eqns  (111)  and  (114) 


x\  _  (  A  —BR~1BT  \  (  x  \ 

x)-[0  -AT  )[x) 


(115) 


with  boundary  conditions  x(t0)  and  A(t;)  =  Sfx(tf).  Writing  out  the  A  part  of  Eqn 
(115)  we  have 

Aj  =  0 

A2  =  0 
A3  =  0 
A4  =  — Ax 

A5  =  —  Aj  (116) 

Ae  =  — A3 
A7  =  tx  X7  —  A4 

Ag  =  TjAg  —  Ag 

Ag  =  Ag  —  Ag 
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Integrating  this  equation  from  tj  to  t  gives 

MO  = 
a2(0  = 
a3(0  = 
a4(0  = 

As(0  = 

a6(0  = 

and  the  resulting  optimal  control  is  from  equation  (79) 

0 

0  O3O3 
*1  t 


S2x2(t/) 

&3X3  (t/) 

-SlXi{tf)(t  -  tf) 
-Wi  {tf){t  -  tf) 
-*3*1  ~  tf) 


u°(t) 

— 

(  0 

0  %T 

1 

can  nov 

K  0  0 

r  be  written 

xx 

= 

x4 

x2 

= 

Xg 

x3 

= 

x6 

X* 

= 

1 . 

x7  —  ~A4  — 

is 

_ 

ri 

xg  -  — A5  = 

Xg 

_ 

r2 

x9  —  —  A6  = 

x7 

= 

r3 

-rxx7 

Xg 

=: 

—r2x3 

X9 

= 

—t2x9 

x{tf) 


r  1 

£2, 

^2 " 

£3, 

r3' 


with  boundary  conditions  at  t0  =  x(0)  and  T  —  tf  —  t.  Integrating  the 
equations  from  t0  to  t  gives 

x7(t)  =  e~T^t~t^x7(t0) 

*•(0  =  e~T,{t-t9)x*{t0) 

x9{t)  =  c-r*(‘-**)x9(f0) 


The  next  set  of  differential  equations  can  then  be  written  as 

x4(0  =  e~Tl(t~t°'>x7(t0)  +  ~Xj {tf) (£  -  t0)  -  {tf) (tf  -  t0) 

x6(0  =  e-T3<‘-‘°>x8(f0)  +  j-x2{tj)(t  -  t0)  -  y-x2(tf)(tf  -  t0 ) 

x,(t)  =  +  Jji3(l/)(f  -  !.)  -  -  1.) 
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(117) 


(118) 


(119) 


bottom  three 

(120) 


(121) 


Integrating  Eqn  (121)  from  t0  to  t  gives 

x4(t)  =  ^(1  ~  t~u(t-to))x7{t0)  +  ^-xi(t/)(t  -  t0)2  -  ^x^tfWf  -  t0)(t  -  t0)  +  z4(t0) 

*6(*)  =  :r(l  ~  e-r3(t"‘o))x8(^)  +  y~x2 (tf){t  -  t0 )2  -  j-x2{tf){tf  -  f0)(t  -  t0 )  +  x6(t0)  (122) 

xe(t)  =  :r(l  -  e-r>(t_to))x9(t0)  +  ^x3(t/)(t  -  t0)2  ~  j-x3[tf)(tf  -  t0)(t  -  t0)  +  x6(f0) 

'3  '3  ~3 

The  results  for  Xi(t),x2(t)  and  x3(t)  can  then  be  obtained  by  integrating  the  above 
three  equations  from  t0  to  t  resulting  in 

*.w  =  -  ^(1  -  +  ^*, (*/)(«  -  o3  - 

£-Xi(tf){tf  -  t„)(t  -  t0)2  +  x4{t0)(t  -  t0)  +  a?x(t0) 

*2(0  =  t~^±x8{t0)  -  ^(1  -  c_ra(‘"t'))x8(<o)  +  g^®2(*/)(*  -  t0)3  - 

^rx2{tf){tf  -  t0)[t  -  t0 )2  +  x5{t0){t  -  t0 )  +  x2[t0)  (123) 

*3(0  =  ^*>(‘0)  -  ^2(1  -  +  £*.(*/)(*  -  g3  - 

■~x3(tf){tf  -  t0){t  -  t0 )2  +  x6(t„)(t  -  t0 )  +  x3{t0) 

The  above  results  can  be  put  into  shorthand  notation  for  t  =  tj  ,  and  T  =  t}  —  t0, 
resulting  in 

(^T3  +  l)*i(t/)  =  xl(tp)  +  Tz4(t0)-^(l-e1)*7(<<») 

(y-T3  +  l)x2(t/)  =  x2{t0)  +  Txb{t0)  -  -^(1  -  e2)z8(«o) 

r2  72 

(fly3  +  i)z3(t/)  =  x3(t0)  +  Tx6(t0)  -  ^j(l  -  e3)x8(t0) 

^-T 2x1(tf)  +  x4(tf)  =  x4(t0)  +  i(l  -  c1)z7(f0)  (124) 

T2xi{tf )  +  xB(t,)  =  x5(t0)  +  ^(1  -  e2)x8(t„) 

4~T2x3{tf)  +  x6(tf)  =  x6{t0)  +  1-{1  -  e3)x9[t0) 

~r3  '3 

x7(tf)  =  e~TlTx7(t0)  =  exx7{t0) 

xs[tf)  =  e~T*Txe(t0)  =  e2x8(t0) 

x9(t/)  =  e~T3Tx9(t0)  =  e3x9(t0) 
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and  in  matrix  notation 


where  e;  = 


[  &r3  + 1 

0 

0 

0 

0 

0 

0 

0 

0  ' 

0 

£T3  +  1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

%T3  +  10 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

fr.T1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

'  1 

0 

0 

T 

0 

0  - 

El 

fj 

0 

0 

• 

0 

1 

0 

0 

T 

0 

0 

-El 

t2 

0 

0 

0 

1 

0 

0 

T 

0 

0 

_Ex 
r • 

0 

0 

0 

1 

0 

0  Ex 

0 

0 

(</)  = 

0 

0 

0 

0 

1 

0 

0 

E2 

0 

x{t0) 

0 

0 

0 

0 

0 

1 

0 

0 

Ez 

0 

0 

0 

0 

0 

0 

h 

0 

0 

0 

0 

0 

0 

0 

0 

0 

e2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

«3 

e~r<r  and  £,•  =  1(1  -  e,).  By  rewriting  the  coefficient  matrix  as 

'll  0  0  000000' 

0  22  0  000000 
0  0  33  000000 

41  0  0  100000 

0  52  0  010000 
0  0  63  001000 

0  0  0000100 
000000010 
0  0  0000001 


(125) 


(126) 


(127) 
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its  inverse  has  been  obtained,  resulting  in 


1 

IT 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

71 
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0 

0 

0 

0 

0 

0 

0 

0 

1 

33 

0 

0 

0 

0 

0 

0 

41 

“TT 

0 

0 

1 

0 

0 

0 

0 

0 

0 

52 

~7l 

0 

0 

1 

0 

0 

0 

0 

0 

0 

63 

33 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

(128) 


Post-multiplying  Eqn  (118)  with  Eqn  (128)  and  the  right-hand  side  of  Eqn  (125)  which 
is  Eqn  (126),  substituting  the  actual  terms  from  the  coefficient  matrix  and  letting  the 
current  time  4 1  =  t0,  results  in  the  continuous  optimal  feedback  navigation  law  with 


Ui°(t)  =  - 

u2°{t)  =  - 
u3°(t)  =  - 


3  sxT 

SiT3  +  3  rx 
3  szT 

szT 3  4-  3r2 
3s3  T 

S3T3  4-  3r3 


xx{t)  - 

x*  (0  - 
xz(t)  - 


Z^T2 
sxT^  +  3rj 


*4  (t)  4- 


3J»i(l  ~tx)T 
r^T*  +  3^)' 


3^2 T2  ,  .  (  3^2(1  —  62)2’ 

T3  *5 ~r  m  2(a  1  O- 


s2r3  +  3r2 
3s3T2 
S3T3  +  3r3 


x6(t)  + 


r22(s22’3  +  3r2) 
3s3(l  —  e3)T 


(129) 


s3T3  +  3r3  v  '  S3T3  +  3r3  6V  '  ’  r32(s3T3  +  3r3)-'9V^ 

If  we  let  r  =  r  =  0  and  s  =  1  in  the  above  guidance  laws,  then  u  —  —  3[yi  -f  If 

we  now  let  V  =  closing  velocity  along  the  line-of-sight  and  a  =  z./T  =  the  line-of-sight 
angle,  then  tq  =  —ZVd  ,  which  is  the  well  known  proportional  navigation  law  with  a 
navigation  ratio  of  3. 

3.5  Estimating  Time-to-go  and  Target  Acceleration 

The  necessity  of  knowledge  of  time-to-go  ( Tgo )  arises  from  the  fact  that  during  the 
development  of  the  guidance  law  the  final  time  was  assumed  fixed  while  it  also  assumed 
complete  control  of  all  three  vehicle  acceleration  components.  In  real  life,  both  hose 
assumptions  will  have  to  be  evaluated  for  application  of  navigation  laws  to  various 
vehicles.  When  the  rocket  engine  thrust  magnitude  is  fixed,  the  solution  of  Tgo  is 
unique.  When  the  rocket  engine  is  throttleable  there  is  a  certain  latitude  in  the  choice 
of  Tgo  for  there  is  then  additional  flexibility  in  the  allocation  of  thrust  acceleration. 
A  Tgo  algorithm  can  then  be  developed,  assuming  zero  target  acceleration,  and  then 
be  expanded  to  include  target  acceleration.  Furthermore,  some  information  such  as 
target  acceleration  for  terminal  navigation  implementation  is  not  directly  available  and 
requires  the  development  of  an  on-board  extended  Kalman  Filter. 

An  important  issue  which  needs  to  addressed  concerns  the  interfacing  of  the  navigation 
and  estimation  algorithms  with  the  terminal  sensors,  the  airframe  sensors,  and  the 
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autopilot.  The  terminal  sensor  measurement  with  the  sensor  gimbal  angles  provide 
information  referenced  to  body  coordinates  and  are  corrupted  by  noise.  The  autopilot 
uses  the  navigation  commands  after  transformation  from  the  inertial-  to  the  body  fixed 
frame  to  maneuver  the  airframe,  which  is  sensed  by  on-board  dynamic  sensors  and  is 
fed  back  to  the  autopilot,  to  provide  system  stability.  This  motion  is  also  sensed  by  the 
terminal  sensor  thus  closing  the  kinematic  loop.  The  dynamic  sensors  for  a  similar  set 
are  also  used  to  update  the  on-board  inertial  reference  frame  used  for  the  navigation 
equation.  Thus  the  coordinate  frame  transformation  between  the  terminal-sensor,  the 
body  frame,  and  the  inertial  frame  need  to  be  considered. 

One  of  the  most  important  considerations  in  the  design  of  homing  missile  systems  is 
the  effect  of  body  motion  coupling  onto  the  signal.  Body  motion  coupling  arises  mainly 
from  the  distortion  of  the  incoming  target  return  energy  as  it  passes  the  radome  which 
protects  the  sensor  and  antenna  assembly  in  the  nose  of  the  missile.  This  issue  is  further 
addressed  in  the  next  section. 

8.6  Radome  Boresight  Errors 

One  of  the  main  concerns  that  must  be  addressed  when  designing  sensor  hardware 
and  guidance  laws  for  homing  missiles  is  that  of  radome/optical  window  boresight  error 
interactions.  Boresight  errors  arise  from  the  distortion  of  the  incoming  energy  wavefront 
as  it  passes  through  the  missile’s  radome/optical  window. 

INCIOENT  \ 

WAVEFRONT  y 


\ 


Fig.  3:  Wave  Mechanism  in  Hollow  Dielectric  Shell 

The  principal  contributor  to  this  distortion  is  the  non-symmetrical  refraction  of  the 
wavefront  that  results  when  the  energy  confronts  the  different  dielectric  properties  at  the 
interface  between  the  atmosphere  and  the  non-hemispherically  shaped  radome.  Typical 
wave  mechanisms  in  a  hollow  dielectric  shell  are  shown  in  Figure  3.  Other  contributors 
to  wavefront  distortion  are  internal  reflections  within  the  radome,  channeling  of  energy 
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along  the  radome  surface,  and  multipath  and  shadowing  caused  by  other  hardware  near 
the  nose  of  the  missile. 


The  ultimate  consequence  of  this  wavefront  distortion  is  an  apparent  tracking  error 
induced  by  the  radome/optical  window.  Hence  the  missile/target  line-of-sight  direction 
which  is  measured  by  the  missile  seeker  appears  to  be  perturbed  from  its  true  value. 
This  perturbation  is  a  strong  function  of  the  direction  with  respect  to  the  radome  from 
which  the  energy  is  coming,  and  thus  is  a  function  of  the  seeker  look  angle  which  in  turn 
is  a  function  of  the  pursuit  geometry  as  it  evolves  over  time.  This  dependence  on  the 
geometry  creates  body  motion  coupling  that  can  alter  the  overall  response  of  the  missile 
system  and  modify  the  measured  noise  characteristics  as  discussed  below. 

The  sensor/missile/target  angular  relationships  are  shown  in  Figure  4,  and  a  simplified 
single-plane  representation  of  a  skid-to-turn  homing  missile  guidance  system  is  shown 
in  Figure  5.  Included  in  this  latter  figure  is  the  body  motion  coupling  feedback  path 
caused  by  radome/optical  window  boresight  interactions.  This  feedback  path  arises 
in  the  missile  as  follows.  The  missile/target  line-of-sight  vector  o  is  tracked  by  the 
missile’s  seeker  and  the  output  of  the  associated  electronics  forms  an  estimate  of  the 
rotation  rate,  at,  of  the  line-of-sight  vector.  This  measured  line-of-sight  rotation  rate 
provides  the  principal  guidance  information  for  the  homing  missile. 


However,  because  of  the  radome  boresight  error  the  measured  line-of-sight  vector  is 
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perturbed  from  its  true  direction  o  by  an  angle  ej.  Thus  the  apparent  line-of-sight 
angle,  ot,  can  be  written  in  terms  of  the  true  line-of-sight  angle,  o,  and  the  radome 
induced  error  as  ot  =  a  +  and  the  apparent  line-of-sight  rate,dv,  is  then: 

dot  _  do  dex  d(0  +  e) 
dt  ~  dt  d(0  +  e)  dt 

or 

d 1  =  0  +  r(/?  +  e)  (131) 

where  r  is  known  as  the  boresight  error  slope,  i.e.  the  change  of  boresight  error  for  a 
change  in  seeker  gimbal  look  angle. 


Fig.  5:  Missile  Guidance  System  Showing  Boresight  Error  Feedback 

But  (3  =  «  =  b  —  xp,  where  d  is  the  look  angular  rate,  and  t/>  is  the  missile  body  rotation 
rate.  Hence  the  measured  line-of-sight  rate  can  thus  be  expressed  as 

dt  =  (1  +  r)<7  +  r\j>  (132) 

This  equation  illustrates  the  major  effects  of  radome/optical  window  boresight  errors 
on  the  measured  guidance  signal.  First  the  radome  error  modifies  the  gain  from  the 
true  guidance  signal,  d,  to  the  measured  guidance  signal,  dt.  Second,  and  by  far  the 
most  important  effect  is  that  it  couples  body  motion,  t />,  onto  the  measured  guidance 
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signal,  at.  The  presence  of  this  terms  results  in  the  closing  of  an  outer  feedback  path 
around  the  airframe,  identical  to  the  standard  rate  gyro  feedback  path  used  in  autopilot 
design  to  achieve  airframe  stability,  and  thus  can  alter  the  system  response  and  stability 
characteristics.  In  particular,  as  the  value  of  r  becomes  more  negative,  the  effective 
missile  system  time  constant  decreases.  On  the  other  hand,  as  r  becomes  more  positive, 
the  missile  system  time  constant  increases,  resulting  in  a  more  sluggish  system  response. 

There  are  a  number  of  loop  parameters  which  affect  the  sensitivity  of  the  system  to  the 
positive  and  negative  values  of  the  boresight  error  slopes.  In  the  guidance  computer, 
both  the  guidance  gain  and  the  guidance  time  constant  affect  the  neutral  stability  slopes. 
If  the  guidance  gain  is  decreased,  the  body  motion  coupling  gain  also  decreases,  and  a 
larger  value  of  |r|  can  be  tolerated.  Increasing  the  guidance  filter  time  constant  decreases 
the  loop  gain  at  high  frequencies.  Hence,  increasing  the  guidance  filter  time  constant 
increases  the  neutrally  stable  boresight  error  slope. 

A  parameter  which  typically  can  not  be  modified  during  guidance  system  design,  but 
which  also  greatly  affects  missile  behavior  is  the  effective  airframe  lift  coefficient  A  (in  g’s 
per  degree  of  angle-of-attack).  In  general,  as  the  magnitude  of  the  effective  lift  coefficient, 
A ,  decreases,  the  body  motion  feedback  loop  gain  associated  with  the  airframe  block  in 
Figure  5  increases.  Consequently,  for  low  lift  airframes  the  effect  of  the  radome  boresight 
errors  on  missile  performance  characteristics  becomes  more  pronounced. 

When  a  radome  induced  instability  occurs  in  a  missile  an  oscillation  begins  to  develop  in 
the  missile’s  steering  control  channels.  This  can  lead  to  a  number  of  phenomena  which 
degrade  missile  performance: 

1.  When  the  boresight  error  slope  is  negative  the  instability  exists  at  a  low  frequency 
well  within  the  guidance  bandwidth  of  the  system.  The  missile’s  response  to  the 
undesirable  steering  command  is  thus  unattenuated  and  the  missile  may  deviate 
radically  from  its  desired  trajectory,  resulting  in  large  miss  distances. 

2.  If  an  instability  occurs  due  to  a  large  positive  boresight  error  slope,  it  results  in 
an  oscillation  which  is  of  high  enough  frequency  that  it  is  typifcally  attenuated  by 
the  missile  airframe  dynamics.  Trajectory  anomalies  are  now  caused  by  dynamic 
range  constraints  within  the  guidance/autopilot  system  and  manifest  themselves 
principally  as  clipping  or  limiting  of  the  oscillatory  guidance  signal.  The  clipping 
results  in  an  effective  reduction  in  the  guidance  gain  from  its  design  value,  thus 
impairing  system  performance. 

An  important  property  of  the  radome  induced  instability  is  that  the  resulting  oscillation 
does  not  necessarily  continue  to  grow  in  amplitude.  Figure  6  illustrates  a  typical  two- 
dimensional  boresight  error  characteristic  as  a  function  of  the  geometric  look  angle,  0. 
Suppose  that  nominally  the  missile/target  line-of-sight  vector  is  offset  from  the  missile 

33 


centerline  by  a  look  angle  of  /30,  and  that  the  magnitude  of  the  boresight  error  slope  r 
at  /?„  is  larger  than  the  critically  stable  slope  rc.  As  a  result,  the  system  will  begin  to 
go  unstable.  The  missile  steering  oscillation  starts  to  grow  and  the  missile  centerline 
begins  to  rotate  relative  to  the  inertially  fixed  missile/target  line-of-sight  vector.  As  the 
look  angle  oscillation  increases  in  amplitude,  averaging  of  the  boresight  error  occurs. 
In  fact,  the  oscillation  will  increase  only  until  the  average  boresight  error  slope  equals 
re.  If  the  oscillation  were  to  increase  further,  the  average  slope  would  be  less  than  the 
critical  slope,  implying  stability,  and  consequently  the  oscillation  would  die  down  until 
the  average  slope  was  again  at  least  re. 


Fig.  6:  Illustration  of  Relationship  between  Boresight  Error 
Characteristics  and  Oscillation  Amplitude 

It  is  through  this  mechanism  that  a  bounded  oscillation  results.  If  the  slope  changes 
are  fast  enough  the  guidance  system  characteristics  will  be  determined  by  an  effective 
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slope  rather  than  the  instantaneous  slope  and  the  two  can  be  substantially  different. 
The  three  parameters  that  need  to  be  considered  are  thus  the  critical  boresight  error 
slope  re,  the  local  or  instantaneous  slope  r,  and  that  local  slope’s  duration  in  look  angle 
A.  In  general,  slopes  that  cover  only  a  small  region  of  the  radome  (i.e.  small  A)  can  be 
tolerated  with  much  larger  magnitudes  than  can  slopes  that  persist  over  a  large  region 
of  the  radome. 

The  presence  of  radome  induced  instability  does  not  automatically  cancel  the  effective¬ 
ness  of  a  missile.  Rather,  if  boresight  error  characteristics  are  such  that  instability 
results  in  bounded  oscillations  of  only  small  amplitudes,  performance  may  still  be  quite 
satisfactory.  Low  frequency  (negative  r  induced)  trajectory  deviations  may  be  small  and 
high  frequency  (positive  r  induced)  guidance  signal  limiting  may  not  be  severe.  ATS 
has  taken  advantage  of  this  fact  and  developed,  using  dual-input  describing  function 
techniques,  a  method  of  boresight  error  compensation  based  on  measured  amplitudes 
and  frequencies. 

3.7  Terminal  Navigation  Law  Evaluation 

Determination  of  the  terminal  navigation  law  in  terms  of  miss-distance  would  seem 
to  require,  in  general,  the  complete  solution  of  the  trajectory  equation  subject  to  the 
appropriate  driving  and  initial  conditions.  For  practical  control  systems,  however,  it  has 
been  found  impossible  to  obtain  such  solutions  analytically.  The  difficulty  arises  from 
the  fact  that  even  a  linearized  trajectory  equation: 

1.  has  time- varying  coefficients, 

2.  may  be  of  quite  a  high  order,  the  order  being  one  higher  than  the  number  of  time 
lags  in  the  control  system. 

Thus  one  has  to  resort  to  simulation  techniques  to  obtain  miss-distance  performance 
results.  Even  then,  because  of  the  time-varying  coefficients,  the  cause-effect  relationship 
is  a  function  of  two  variables  (one  being  the  time  of  application,  r,  of  the  cause  and  the 
other  being  the  time  of  observation,  t  of  the  effect)  and  a  straight  forward  simulation  has 
to  be  run  many  times,  each  for  a  different  maneuver  application  time,  r,  which  is  a  very 
costly  and  time-consuming  process.  American  Technical  Services  (ATS)  has  overcome 
the  computational  difficulty  with  the  use  of  an  adjoint  miss-distance  simulation.  The 
state  transition  matrix  of  an  adjoint  system  is  related  to  the  state  transition  matrix  of 
the  original  system  by  an  interchange  of  the  running  time,  t,  and  application  time,r,  and 
a  transposition.  Thus,  while  the  behavior  of  the  system  with  respect  to  t  is  a  function 
of  the  dynamics  of  the  original  system,  the  behavior  of  the  system  with  respect  to  r  is 
a  function  of  the  dynamics  of  the  adjoint  system.  As  a  result,  the  adjoint  simulation 
requires  only  a  single  run  to  produce  terminal  miss-distance  results  due  to  each  of  the 
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initial  conditions  and  each  of  the  forcing  functions  and  system  disturbances.  The  adjoint 
technique  of  analysis  is  further  discussed  in  Appendix  A.  In  that  appendix,  Figure  A-l 
shows  a  block  diagram  of  ATS’s  forward  simulation  with  initial  turning  rate,  heading 
and  pointing  errors,  as  well  as  target  maneuver  inputs  and  noise  disturbances.  Figure 
A-2  shows  the  block  diagram  of  the  corresponding  adjoint  simulation  with  outputs  in 
terms  of  miss-distance  due  to  initial  turning  rate,  heading  and  pointing  errors,  as  well 
as  miss-distances  due  to  each  of  the  maneuvers  and  noise  disturbances. 

The  terminal  navigation  law  evaluation  was  performed  using  the  adjoint  simulation.  The 
adjoint  method  of  analysis  is  especially  valuable  when  the  system  is  subject  to  multiple 
input  disturbances,  since,  as  discussed  in  Appendix  A,  the  sensitivities  of  the  system 
output  (i.e.  terminal  miss-distance)  to  the  specific  input  disturbances  are  obtained  with 
a  single  run  at  the  outputs  of  the  adjoint  system. 


Noise,  and  target  glint  noise  in  particular,  has  long  been  one  of  the  primary  concerns  in 
assessing  guidance  performance.  Noise  interactions  with  ra dome  boresight  errors  have 
not  received  much  attention,  since  with  a  forward  simulation  it  requires  many  runs, 
often  hundreds,  to  obtain  statistically  meaningful  results.  The  semi-active  receiver  noise 
model  used  in  the  adjoint  simulation  was  based  on  the  rms  angular  noise  equation 


knJzsjif 


(133) 


and  on  the  fact  that  the  noise  power  density  (0(o))  times  the  bandwidth  equals  the 
angular  error  squared  (ae2  —  6 (o)  *  Bu),  then 
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The  form  of  the  radar  equation  used  for  computation  of  the  signal  noise  ratio  is 

C  /  J\T  _  _ PgGjGyjX^OT _ 

'  ~  ( 4n)*kTNTLRIT2RMT2Bu  {  } 

and  thus 

<(.)  =  (136) 

2km2  PaGTGw\2oT  '  " 

The  range  independent  noise  model  used  in  the  adjoint  simulation  used  a  noise  power 
density  of 

«(»)  =  ^(2J5rn)!Ar  rad’/Hz  (137) 

The  glint  noise  model  used  is  based  on  a  model  for  a  target  that  consists  of  scatterers 
distributed  uniformly  over  the  target  length  and  is  given  by 


a gi  =  0.35L  feet 


(138) 


Target  maneuver  models  consist  of  a  target  step  displacement,  a  target  step  velocity 
and  a  target  step  acceleration,  a  filtered  step  acceleration,  and  a  sinusoidal  target  accel¬ 
eration.  The  initial  condition  errors  consist  of  initial  heading  error  turning  rate,  initial 
missile  turning  rate,  and  seeker  pointing  error.  The  miss  distance  results  for  zero  bore- 
sight  error  slope  due  to  target  maneuver,  initial  heading  error  and  initial  turning  rate 
are  shown  in  Figure  7. 

A  boresight  error  model  based  on  Eqn  (132)  was  then  implemented  in  the  adjoint  sim¬ 
ulation  and  the  miss  distance  sensitivity  boresight  error  slopes  determined.  They  were 
obtained  by  keeping  the  noise  inputs  constant  and  increasing  the  boresight  error  slopes 
for  an  rms  miss  distance  due  to  glint  noise  twice  the  original  value.  The  miss  distance 
sensitivities  of  the  glint,  amplitude  and  receiver  noise  are  shown  in  Figure  8. 

The  parameter  values  used  in  Eqns  (134-139)  for  the  missile  are  0B  =  .72  deg,  kT  ~  4 
x  10~21  w/Hz,  Nf  =  12  db,  L  =  15  db,  km  =  1.25,  PA  =  3  kW  (100  kW  peak),  GT  =  75 
db,  Gw  =  49  db,  A  =  .0015  m,  N  =  15  bits,  R  —  50  degrees  and  A  =  ^j. 

The  rms  miss  distance  result  due  to  the  semi-active  receiver  noise  for  a  missile-target 
closing  velocity  of  10  kft/sec,  a  time  constant  of  .1  sec  and  an  illuminator-target  range  of 
400  mi,  is  shown  to  be  at  .45  feet.  The  rms  distance  results  due  to  the  range  independent 
noise  is  .48  feet. 

For  the  major  length  (32.8  ft)  of  the  target,  the  resulting  rms  miss  distance  due  to  glint 
noise  was  8.81  ft  and  for  the  minor  axis  (9.84  ft)  of  the  target,  the  resulting  rms  miss 
distance  was  found  to  be  1.64  feet. 

The  total  rms  miss-distance  due  to  all  three  noise  sources  (i.e.  receiver,  range  indepen¬ 
dent  and  glint  noise)  can  be  obtained  by  obtaining  the  root-sum-squared  value  and  is 
rms  miss  =((.45)2  +  (.48)2  -I-  (8.81)2)a  =  8.83  ft  for  the  major  target  dimension  and  rms 
miss  =  ((.45)2  +  (.48)2  +  (2.64)2)?  =  2.72  ft  for  the  minor  target  dimension.  Due  to  the 
ratios  of  the  glint  miss  to  receiver  or  independent  miss,  the  resultant  total  rms  miss  is 
essentially  the  miss  distance  due  to  the  glint  noise. 

Not  taken  into  account  in  the  above  noise  models  has  been  the  degree  of  correlation  of 
the  error  or  the  amplitude  from  one  observation  to  the  next.  In  the  actual  missile  this 
can  be  important  because  it  will  determine  the  response  of  the  measurement  system  as 
well  as  the  effectiveness  of  a  smoothing  filter.  Decorrelation  may  be  caused  by  target 
rotation  relative  to  the  line-of-sight  or  intentionally  created  using  a  frequency  diversity 
system. 

To  achieve  a  hit  capability  of  more  than  or  equal  to  99.5%,  it  is  required  that  o2mi„  < 
(£)2,  where  L  is  the  target  dimension.  Thus  for  a  hit  capability  equal  to  99.5%,  the  rms 
miss  for  the  target  dimension  is  (^)  =  5.47  ft  and  for  the  minor  target  dimension  is 
(S^4)  =  1.64  feet 
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Since  the  present  miss  distance  for  the  weapon  is  essentially  due  to  glint  noise,  it  can  be 
reduced  to  the  desired  values  for  a  PH  of  99.5%  by  either  the  use  of  frequency  diversity 
or  by  increasing  the  weapon  response  time.  Although  increasing  the  weapon  response 
time  will  decrease  the  miss  due  to  glint  noise,  it  will  increase  the  miss  distance  due  to  all 
other  sources.  The  use  of  frequency  diversity  is  therefore  considered  to  be  the  preferred 
method  of  decreasing  the  miss  distance  and  has  been  analyzed  below. 

The  number  of  independent  samples  required  to  achieve  a  PH  =  99.5%  is 

__  .miss  achieved.2  _  ,2.64. 2 ,8. 81, 2  _ 

^  —  '  miss  desired  '  '1.64'  '5.47'  ~ 

For  uniformly  distributed  scatterers,  the  correlation  frequency  interval  fe  is  given  by 

/‘  =  2i  =  Tr^r  =  50  MH2  («°> 

where  c  is  the  speed  of  light  and  L  the  minor  target  dimension  in  meters.  The  required 
bandwidth  A  /  can  be  determined  from 

A/  =  (rif  -  l)fe  =  1.59  *  50  =  80  MHz 

which  is  .04%  of  the  200  GHz  radar. 

Miss  distances  versus  time-to-go,  due  to  transient  type  causes  are  shown  in  Fig.  7  for  an 
initial  heading  error  of  one  degree,  a  filtered  step  acceleration  of  32.2  feet/sec2  through  a 
2  second  time  constant,  a  sinusoidal  target  acceleration  of  32.2feet/sec2  at  a  frequency  of 
—5  rad/sec  and  of  a  step  target  displacement  of  100  feet.  It  is  clear  from  the  figure  that 
a  step  target  acceleration  just  before  intercept  or  at  times-to-go  of  .65  and  1.7  seconds 
will  result  in  the  maximum  benefit  to  the  target.  On  the  other  hand,  any  maneuvers 
before  a  time-to-go  of  three  seconds  and  at  .28  and  1.15  seconds  to  go  are  of  no  benefit 
to  the  target. 

3.8  Navigation  Law  Implementation 

Implementation  of  the  navigation  laws  will  be  relatively  simple  using  the  on-board  avail¬ 
able  sensors  plus  an  extended  Kalman  filter  for  the  terminal  navigation  law,  to  estimate 
target  motion  and  time-to-go.  Because  the  measurements  from  the  sensor  are  in  the 
sensor’s  coordinate  system,  the  body  motion  sensors  are  in  the  vehicle’s  body  fixed  co¬ 
ordinate  system,  and  the  navigation  commands  to  the  autopilot  have  to  be  in  the  vehi¬ 
cle’s  body  fixed  coordinate  system,  it  is  necessary  to  implement  transformation  matrices 
relating  the  sensor’s  coordinate  system  to  the  vehicle’s  body  fixed  coordinate  system  to 
the  inertial  coordinate  system.  Since  all  those  equations,  including  the  Kalman  filter 
and  navigation  laws,  are  quite  readily  expressed  in  terms  of  vectors  and  matrices,  and 
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since  their  calculations  are  naturally  performed  by  matrix  manipulations,  it  becomes 
rather  attractive  to  provide  the  navigation  computer  with  an  interpreter  that  translates 
a  powerful  convenient  set  of  matrix  and  vector  instructions  into  machine  language  at 
the  time  of  machine  execution  of  each  instruction.  If  all  those  computations  were  pro¬ 
grammed  as  scalar  equations,  a  great  deal  of  storage  for  programs  would  be  required 
and  a  great  deal  of  the  elegance  of  the  equations  would  be  lost.  On  the  other  hand,  an 
interpreter  requires  some  storage  and  could  increase  the  execution  time  of  the  navigation 
computations,  which  in  turn  can  affect  the  accuracy  and  stability  of  the  navigation  and 
control  loop. 

When  the  rocket  engine  is  throttleable,  there  is  a  certain  amount  of  latitude  in  the  choice 
oiTgo  because  there  is  then  additional  flexibility  in  the  allocation  of  thrust  acceleration. 
The  navigation  solutions  determine  the  required  allocation  of  thrust  acceleration  along 
each  coordinate  axis.  Since  the  component  of  thrust  acceleration  along  any  axis  is 
obtained  from  a  single  thrust  acceleration  vector,  the  total  thrust  acceleration  must 
be  time-shared  among  the  controlled  axes.  The  program  for  time-sharing  this  thrust 
acceleration  must  be  carefully  planned  so  that  the  solutions  to  the  problems  along  the 
controlled  axes  are  realized  simultaneously. 
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Fig.  8:  Miss  due  to  Noise  Causes  for  a  Sampling  Rate  of  40  cps 


4.0  Conclusions 


The  navigation  laws  developed  in  this  report  can  control  final  coordinates  of  positions 
as  well  as  final  components  of  velocity  and  control  throttleable  as  well  as  fixed  thrust 
rockets.  The  steering  commands  are  directly  in  terms  of  the  current  and  desired  bound¬ 
ary  values  of  the  components  of  the  position  and  velocity  vectors.  The  equations  are 
exact  and  can  accommodate  any  gravitational  field  model.  By  treating  the  weighting 
parameters  explicitly,  universal  navigation  laws  have  been  obtained  that  are  applicable 
to  many-faceted  complex  missions. 

For  applications  to  the  various  missions,  the  weighting  parameters  should  be  optimized 
and  unified  as  much  as  possible,  taking  into  consideration  the  actual  fuel  budgets,  fixed 
or  throttlable  rockets,  vehicle  dynamics  and  other  factors,  to  establish  the  optimal  levels 
of  control  for  those  missions.  The  most  economical  and  effective  method  to  evaluate  the 
algorithms  for  any  mission  is  to  implement  a  detailed  simulation  of  the  system  and  to 
perform  simulated  fly-outs  against  maneuvering  and  non-maneuvering  targets. 

Such  a  simulation  could  consist  of  a  six-degree-of-freedom  model  of  the  system,  contain¬ 
ing  detailed  math  models  of  the  major  subsystems,  including  the  sensor,  autopilot  and 
propulsion,  detailed  aerodynamic  models  of  the  airframe  characteristics  supported  by 
wind  tunnel  generated  aero  data,  and  the  models  that  describe  the  vehicle’s  equation 
of  motion.  In  addition,  the  simulation  could  contain  a  three-degree-of-freedom  target 
model  that  incorporates  a  suitable  evasive  maneuver  algorithm.  To  obtain  actual  sys¬ 
tem  performance  the  transfer  alignment  errors  at  initialization  and  the  instrument  drifts 
during  flight  need  to  be  implemented.  They  will  affect  the  size  of  the  error  basket  at 
handover  to  terminal  navigation.  Once  in  the  terminal  navigation  mode,  those  errors 
are  of  no  further  concern. 

For  the  terminal  navigation  law  using  a  terminal  sensor,  an  extended  Kalman  filter, 
which  is  typically  based  upon  the  same  equations  of  motion  as  used  for  the  terminal 
navigation  law,  should  be  developed  and  implemented  in  the  simulation.  In  addition,  this 
filter  should  also  be  used  to  obtain  an  estimate  of  time-to-go.  During  the  development 
of  the  Kalman  filter  it  will  also  be  necessary  to  investigate  in  which  coordinate  reference 
frame  (sensor,  body,  or  inertial)  the  Kalman  filter  should  be  implemented. 

Mechanization  of  the  navigation  laws  and  Kalman  filter,  together  with  the  coordinate 
reference  transformations,  should  then  be  implemented  and  evaluated  on  the  simula¬ 
tion  which  will  lead  to  the  navigation/estimation  mechanization  specifications  for  the 
navigation  computer. 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


H  -  Hamiltonian 

J  -  performance  criterion/optimal  control  trajectory 

L  -  minor  target  dimension  in  meters 

M  -  missile  position  vector 

T  -  time  to  orbit 

a  -  acceleration 

c  -  speed  of  light 

g  -  gravitational  force 

h  -  constant  specific  angular  momentum  of  the  orbit 

r  -  position  vector  from  the  center  of  the  earth 

u  -  control  vector  of  acceleration 

v  -  velocity  vector 

z  -  state/costate  vector 

Ui°(t)  -  continuous  optimal  feedback  control  terms 

$  -  transition  matrix 

^  -  angle  between  r  and  an  arbitrary  reference 

6f  -  bandwidth 

A  -  costate  vector  of  the  Lagrange  multiplier 

/j.  -  constant,  defining  the  specific  gravitational  force 

t  -  true  line-of-sight  angle 

Tt  -  apparent  line-of-sight  angle 

Tl  -  rotation  rate  of  apparent  line-of-sight  angle 

IRU  -  Inertial  Reference  Unit 

BIM  -  Ballistic  Intercept  Missile 

INS  -  Inertial  Navigation  System 

Tgo  -  Time-to-go 

FOV  -  Field  of  view 
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APPENDIX  A 

ADJOINT  METHOD  OF  ANALYSIS 


The  adjoint  method  of  analysis  has  proven  to  be  a  very  use¬ 
ful  concept  for  studying  the  effects  of  forcing  functions  and 
initial  conditions  on  linear  combinations  of  the  state  variables 
of  the  original  system. 

If  the  original  time-varying  system  can  be  described  by  the 
vector-matrix  equation: 

x(t)=A(t)x(t)  A-l 

then  the  corresponding  adjoint  system  is  defined  by 

a  ( t )  =  -A  (  t)Ta(t)  A- 2 

The  solution  of  the  original  system  is  given  by 
X(t)*$>(t,T)X(T)  A- 3 

where  <f  (t,x)  is  the  state  transition  matrix  and  satisfies 

$(t,T)  =  A(  t )  <M  t ,  t  )  ,  $  (  x  r  x )  =  I  A-4. 

Similarly,  the  solution  to  the  adjoint  system  is  given  by 

a(t)*\jj(t,T)a(T)  A- 5 

where  t(t,T)  is  the  state  transition  matrix  and  satisfies 

'J(t,T)  =  -AT(t)\p(t,T)  ,  tp(  t,  t  )  =  I  A-6 

The  state  transition  matrix  for  a  time-varying  system  has 
two  variables,  one  being  the  application  of  the  cause,  and  the 


other  the  time  of  observation  of  the  effect.  The  usefulness  of 
the  adjoint  system  stems  from  the  fact  that 


*T<t,x)  =  C  $ 

A- 7 

^T(t,x)  -  4>(  x,t) 

A-8 

Thus  the  behavior  of  the  system  with  respect  to  the  variable  t  is 
a  function  of  the  dynamics  of  the  original  system  while  the 
behavior  of  the  system  with  respect  to  the  second  variable  t  is  a 
function  of  the  dynamics  of  the  adjoint  system. 

The  vector-matrix  equation  for  a  system  with  input  u(t)  and 
output  y(t)  is  given  by 

x(t)  =  A(t)xit)  +  B(t)u(t)  A-9 

ytt)  =  C(t)x(t)  A- 10 

where  A(t)  is  an  n  x  n  system  matrix 
B(t)  is  an  n  x  m  input  matrix 
C(t)  is  a  p  x  n  output  matrix 
and  the  system  has  m  inputs  u(t)  and  p  outputs  y(t). 

The  general  solution  to  the  above  system  is  given  by 

x(t)  =  $  ( t  ,t  )  x(x  )  +  /V  t ,  v)  B(  v )  u(  v)dv  A-ll 

'X 

t 

y(t)  *  C(  t)$  ( t,x  )x(x  )  +/  C(  t  )<{>  ( t ,  v)  B(  v)u(  v)dv  A-12 

where  4>(t,  x)  is  as  before  the  state  transition  matrix  of  the 
homogeneous  equation  x(t)  =  A(t)  x(t)  and  satisfies 

$(t,x)  =  A(  t  )3>  ( t  ,x  )  ,  =  I  A-13 

The  adjoint  system  corresponding  to  the  original  time-varying 
system  is  again  given  by 

a  ( t )  =  -  A  Ta(  t )  A-14 


whose  solution  is  given  by 
a(t)  »  <J/(t,x)  a  (t  ) 


A-15 
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where>Mt,T)  satisfies  as  before 

<I»(t,T)  *  -AT  (t)  i|»(t,T>  A-16 

Taking  the  inner  product  of  the  original  system  solution  (A-ll) 
with  the  adjoint  solution  (A-15)  gives 

t 

cJ(  t )  x(  t )  =  aT(x  )^T(t,T)<£(  t  ,t  )  x(  x  )  +  v)^T(  t,v)$(t,v)B(v)u(v)dv  A- 17 

T 

or  after  using  A- 8 

c^ttJxCt)  =c?<t)x(t)  +  v)B(  v)u(  v)dv  A-18 

For  a  fixed  terminal  time  T  and  for  t  =  -  » 

T 

<f(T)  x(T)  B( v )  u(v)dv  A-19 

since  a ( -« )  =  0 . 

Equation  A-19  can  now  be  used  to  determine  the  linear  combi¬ 
nations  of  the  state  variables,  xi(t)'s,  once  the  adjoint  equa¬ 
tion  has  been  solved  for  the  a(v)'s. 

In  the  time-varying  case  this  generally  must  be  done  by 
simulation  because  of  the  difficulty  in  analytically  solving 
time-variable  differential  equations. 

Solution  of  A-19  requires  specification  of  the  boundary  con¬ 
ditions  on  a(v).  The  boundary  conditions  which  should  be  speci¬ 
fied  depend  upon  the  problem  to  be  solved.  For  example,  if  the 
effect  of  the  forcing  functions  and/or  initial  conditions  on  x1 
(t)  is  to  be  determined,  the  appropriate  boundary  conditions  are 

aj(t)  =6  and  a2^)  =a^(t)  =  . an(t)  =  0. 

where  6  is  the  Kronecker  delta  or  unit  impulse.  Thus  a 

judicious  selection  of  the  terminal  constraint  can  easily 

result  in  A-19  containing  oq( T ) X]( T )  as  the  only  component  at  time  T. 
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If  a(T)  is  chosen  as 

aT(T)  =  CqL(T)  C^T)  Ci3(T)  . Cin(T>:i  A-20 

where  the  C  i  ' s  are  the  ith  row  elements  of  the  p  x  n  output 

matrix  C  at  terminal  time  T,  then  the  output 

y,(T)  =  CC..(  T)  C.0(  T )  . C,  (T)3  x(T)  A-21 

1  i  ll  i2  in 

can  be  obtained,  or 
T 

y.  (T)  =  f  aT(v)B(v)u(v)dv  A-22 

In  order  to  make  these  boundary  conditions  initial  conditions  in 
the  simulation  let  v  =  T  -x,  or  x,=  T  -  v. 


Since  v  runs  "forward  in  time”  from  -  °°to  T,  xqruns  "backward 
in  time"  from  x  *  0  to  Xj*  00  ,  thus 


y±(T)  =  -  ^a^T-x^BCT-T^utT-TpdT!  A-23 

and  the  corresponding  adjoint  system  is 

^(T-Vad-Tj)  A.24 

The  coefficient  matrix  of  this  equation  is  the  transposed  A 
matrix  of  the  original  system,  with  t  replaced  by  T  -tj.  The 
initial  conditions  for  this  system  are  then 

aT  (  0 )  =  CCijfT),  Ci2<T)  . ...  Cfa(  T)D  A-25 

Those  initial  conditions  can  be  established  by  applying  a  unit 
impulse  at  x^  0  through  gains  of 


C 

C 

C 
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(T 

(T 
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) 

) 
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(T  -  Xj) 
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Since  B  and  C  are,  respectively,  the  input  and  output 
matrices  of  the  original  system,  the  net  result  is  that  the  ith 
row  of 

00 

)B(T  -t1)u(T  -TjJdTj  A-27 

is  generated  for  variable  t  by  reversing  the  input  and  output  of 
each  of  the  simulation  elements  of 

x  =  A( t ) x  +  B ( t ) v  A-28 

y  =  C( t )  x 

and  replacing  the  time  t  of  any  time-varying  gains  by  T  -  Tj 
If  the  original  system  has  m  inputs  and  p  outputs,  the  adjoint 
system  has  p  inputs  and  m  outputs.  Each  simulation  run  of  the 
adjoint  simulation  produces  m  outputs.  These  m  outputs  comprise 
a  row  of  C  where  the  running  variable  of  simulation  time  is  T . 
The  particular  row  of  the  matrix  is  determined  by  the  input  on 
which  the  unit  impulse  is  placed  at  the  beginning  of  the  simula¬ 
tion  run. 

For  example,  if  the  unit  impulse  is  placed  on  the  state 
variable  that  represents  miss  distance  in  the  forward  simulation, 
the  outputs  of  the  adjoint  simulation  represent  the  miss  distance 
caused  by  the  corresponding  input  variable  of  the  forward  simula¬ 
tion. 


Starting  with  a  block  diagram  representation 'of  the  forward 
system,  construction  of  the  adjoint  simulation  takes  place  by 

1.  Reversing  all  signal  flows 

2.  Redefining  branch  points  as  summing  blocks 

3.  Redefining  summing  blocks  as  branch  points 

4.  Replacing  t  by  T-  t^in  the  arguments  of  all  time-varying 
coefficients . 
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APPENDIX  B 
ADJOINT  SIMULATION 


In  order  to  generate  the  adjoint  simulation  it  was  first 
necessary  to  develop  a  block  diagram  of  the  forward  simulation 
after  which  the  adjoint  simulation  could  be  developed  using  the 
rules  stated  in  Appendix  A,  i.e. 

1.  Reverse  all  signal  flows 

2.  Redefine  branch  points  as  summing  blocks 

3.  Redefine  summing  blocks  as  branch  points 

4.  Replace  t  by  T  -  in  the  arguments  of  all  time-varying 
coefficients . 

The  space  diagram  showing  the  geometric  relationship  between 
the  missile  and  target  is  shown  in  Figure  B-l.  The  missile  and 
target  velocity,  Vm  and  Vt  respectively,  are  assumed  to  be  con¬ 
stant 

Gravity  effects  are  neglected  and  the  encounter  is  assumed  to  be 
restricted  to  the  x-y  plane.  The  angles  are,  of  course,  subject 
to  change  because  both  missile  and  target  are  free  to  maneuver  in 
the  x-y  plane. 

The  fundamental  relations  governing  the  missile  and  target 
paths  during  the  terminal  engagement  are  the  following  velocity 
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Fig.  B-l :  Space  Diagram  for  a  Homing  Missile 


equations  5 


y  t  =  Vtsin  Y T 

ym  =  VmSinYM 
X  t  =  Vtcos  Y  M 

x  t  =  Vtcos  Y  T 


B-l 

B-2 

B-3 

B-4 


The  angles  YMand  YT  start  at  t  =  0  from  some  initial 


values 


instantaneous  angles  of 


Ytq  and  ymo  and  are  perturbed  by  small  amounts  Ym  and  Yt  during  the 
encounter.  Under  these  conditions ,  the 
the  velocity  vectors  are 

YM  =  YM0  +  Ymand  ^ 

so  that  the  linear  velocity  components  are 


=  yto+  Yt 


yt  =  Vtsin(YT0+Yt)  =  VtsinYT0+YtVtcosYT0 

ym  =  VmSin(YMO+Ym)  =  VinYM0+YmVmCOsYM0 
xt  =  Vtcos(YT0+Yt>  =  VtcosYT0-YtVtsinYT0 

X  =  VmCO8(YM0+Ym)  =VmCOSYM0-YmVmSinYM0 


B-5 

B-6 

B-7 

B-8 

B-9 


m  m  '  MU  m  ra  nu  ra  m 

where  the  small  angle  approximations  sin  y  =  Y  and  cos  Y  3  1  have 


been  used. 
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If  we  assume  that  the  missile  and  target  are  initially  on  a 
lead  collision  course,  y  *  0  and 

VtsinyT0*  Vm  sinYMo  B-10 

and  the  projected  miss  rate  My  is  given  by 

My  »  yt-  ym  KYtVtcosYM0“YmVmcosYMO  B-ll 


Using  (B-8),  (B-9)  and  (B-10)  we  have 

xt-  xm  =  (VtcosYT0-  VmcosYM0>  -  <Yt-Ym)VinsinYM0  B-12 
If  we  neglect  the  difference  term  in  B-12  involving  Yt*  Ym  ,  the 
closing  rate  Vc  is  given  approximately  by 

-  Vc  *  Xt-  Xm  «  V tCOS Yxo"  VmcosYMO  B-13 

and  the  missile-to-target  range  R(t)  is 


The  target 
where 


R(t)  =  R  -V  _t 
o  c 

bearing  as  viewed  by  the  missile  homing 


M 

a  -X 


R(t) 


i_  ftM 
(t)  i  3 


dt. 


B-14 

system  is  cr. 


B-15 


A  block  diagram  representative  of  the  linearized  kinematics  for 
a  homing  missile  is  shown  in  Figure  B-2. 


The  kinematics  represent  that  portion  of  the  problem  that 
one  must  live  with,  i.e.  those  are  the  elements  of  the  problem 
about  which  one  has  no  choice,  and  which  must  be  accepted  the  way 
they  are. 

The  system  must  then  operate  on  the  sight  angle  a  in  such  a 
way  that  the  final  value  of  M  y  (when  t  =R0/vcor  xt  *  xm  )  is 
made  as  small  as  possible. 


Two  observations  can  be  made  at  this  point.  First,  a  double 
integration  appears  in  that  portion  of  the  kinematic  loop  of 


Fig.  B-2;  Linearized  Homing  Missile  Kinematics 

Figure  B-2  joining  ym  to  a.  Consequently,  the  control  portion  of 
the  loop,  whatever  it  may  be,  must  necessarily  provide  some 
derivative  effect  if  system  stability  is  to  be  insured. 
Secondly,  since  the  kinematic  loop  possesses  a  time  dependent 
division  by  range  R(t),  a  linear  time-varying  system  results  even 
ifymand  a  are  coupled  by  a  constant  coefficient  control. 

A  practical  method  of  control  that  has  received  widespread 
application,  and  which  is  used  in  the  HAWK  missile,  is  propor¬ 
tional  navigation.  In  this  type  of  control  one  tries  to  make  the 
missile-path  turning  rate  proportional  to  the  line-of -sight  rate. 
It  is  clear  that  if  the  angular  rotation  of  the  line-of-sight 
between  missile  and  target  can  be  made  to  vanish  and  the  range 
rate  is  of  such  a  sign  as  to  result  in  closure,  collision  is 
inevitable.  Of  course,  the  line-of-sight  rate  cannot  be  reduced 
to  zero  and  held  there  in  any  practical  system  because  of  a 
number  of  system  imperfections,  which  lead  to  uncertainty  in  the 
true  target  bearing. 
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The  true  miss  distance  MT  can  be  defined  as 

Mt  -  Min  /x  t  -  ^  )2  +  t  -  ym  )2  '  B-16 

which  is  essentially  the  distance  of  closest  approach  between  the 
missile  and  target.  The  radical  in  B-16  is  also  the  instantane¬ 
ous  range  between  the  missile  and  target.  In  our  linearized  ver¬ 
sion  of  the  missile  system  we  have  chosen  the  coordinate  system 
in  such  a  way  that  the  approximation  equations  give 


and  used  instead  the  projected  miss  distance  M  defined  by 

My  =  y  t  ‘  yr  0  <,  t  <  R0/Vc  B- 18 

The  forward  missile  simulation  is  shown  in  Figure  B-3.  The  sys¬ 
tem  operates  on  the  sight  line  angle c  ,  derived  from  My  through 
division  by  range  XMT  =  (xt  -  xm>,  in  such  a  way  that  the  final 
value  of  My  (when  t  =  R0/Vcorxt  =  xm)  is  made  as  small  as 
possible. 

The  HAWK  navigation  system  is  implemented  with  a  space- 
stabilized  radar  seeker  aboard  the  missile,  whose  antenna  pro¬ 
vides  a  measure  of  the  target  bearing  a  .  Referring  to  Figure 
B-l ,  the  centerline  of  the  seeker  antenna  makes  an  angle  0  with 
the  reference.  The  angular  error  in  target  bearing  is 

e  -  a  -  0  B-19 

The  angular  error  is  converted  to  a  voltage  by  the  radar  and  fil¬ 
tered  to  provide  the  necessary  driving  signal  to  process  the 
seeker  antenna  in  such  a  direction  as  to  try  continually  to 
reduce  the  angular  error  to  zero.  The  filtered  drive  signal  is 

accordingly 

„  ke 

E  =  -t - ; - 

1  +  IS 
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where  a  simple,  time-constant  lag-filter  has  been  used.  The  pre¬ 
cession  rate  of  the  seeker  9  is  proportional  to  E  ,  so  that 

9  =  kE  B-21 

where  k  is  the  precession  constant  for  the  seeker. 

The  same  drive  signal  E  is  used  as  input  to  the  guidance 
computer  where  the  proportional  guidance  law  and  a  first  order 
noise  filter  are  implemented.  The  output  of  the  guidance  com¬ 
puter  is  the  normal  acceleration  command  to  the  autopilot  which 
can  be  implemented  in  the  simulation  as  a  first,  second  or  third 
order  system. 

The  output  of  the  autopilot  is  the  missile's  normal 
acceleration  in  g's  which  upon  multiplication  by  32.2/Vm  results 

in  the  missile  velocity  vector  rotation  rate  ym  in  rad/sec. 

*  # 

Angle  of  attack  a  is  proportional  to  ym,  i.e.  a  =  kaYm,  where  k 
is  the  inverse  of  the  effective  lift  coefficient  A.  The  control 
portion  of  the  missile  homing  system  is  completed  by  integrating 
the  missile  velocity  vector  rotation  rate.  Thus  ym  together  with 
possible  target  maneuvers  are  input  to  the  homing  kinematics  loop 
whose  output  after  division  by  XMT  forms  the  input  to  the  forward 
simulation. 

Boresight  error  slope  is  implemented  by  forming  a-\p  and  sub¬ 
tracting  that  signal  from  a  after  mult iplicatioa  by  kye  =  -r  , 
thus  creating 

a'  *  a  +  r(o-<p)  =  (1  +  r)a  -  r ip  B-22 

The  output  quantity  of  interest  is  My,  the  projected  miss 
distance.  XMT  is  the  range  at  any  time  t  after  the  start  of  the 
homing. 


The  initial  geometry  is  fixed  by  the  missile  and  target  heading 
at  the  start  of  the  homing. 

The  glint  or  scintillation  noise  input  is  represented  as  a 
linear  displacement  at  the  target  and,  like  miss  distance,  is 
divided  by  range  to  become  an  angular  noise  input  to  the  seeker. 
To  allow  target  glint  representation  as  stationary  band-limited 
noise  with  amplitude  and  bandwidth  dictated  by  the  physical 
dimensions  and  motion  of  the  target,  a  band-limited  glint  noise 
input  is  provided.  Since  the  statistical  properties  of  band- 
limited  noise  is  the  same  as  filtered  "white"  noise,  a  simple  lag 
filter  between  the  band-limited  input  and  the  driving  white  noise 
is  required.  The  amplitude  or  fading  noise  is  introduced  as  an 
angular  uncertainty  one'  and  is  assumed  to  be  white.  Semi¬ 
active  receiver  noise  is  introduced  at  the  same  place  after  mul¬ 
tiplication  by  XMT  squared. 

Consider  now  the  initial  condition  inputs  to  the  system 
indicated  by  double  arrows  in  Figure  B-3.  One  such  input  is  the 
head  pointing  error  at  the  start  of  homing  and  is  assumed  to  be  a 
constant  random  variable.  The  initial  turning  rate  and  heading 
error  of  the  missile  is  shown  as  an  input  to  Ym  and  Ym  respec¬ 
tively.  Either  of  these  errors  could  be  zero  under  certain  cir¬ 
cumstances.  If  the  seeker  is  locked-on  and  settled  out  at  start 
of  homing,  the  seeker  initial  error  input  becomes  zero.  If  the 
missile  is  not  maneuvering  at  the  start  of  homing,  the  initial 
turning  rate  and  heading  error  al30  vanish.  Target  maneuvering 
inputs  consisting  of  step,  filtered  step  and  sinusoidal  target 
maneuvers  are  also  sown  on  Figure  B-3. 

The  forward  simulation  has  12  inputs  and  a  single  output, 
miss  distance.  The  system  adjoint  must,  therefore,  have  a  single 
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input  and  12  outputs.  The  adjoint  is  generated  by  reversing  all 
inputs  and  outputs  and  generating  all  time  functions  backwards. 
The  adjoint  missile  simulation  is  shown  in  Figure  B-4.  It  has  a 
single  impulse  input,  indicated  by  the  double  arrow  in  Figure  B- 
4,  at  the  point  where  the  miss  distance  was  measured  in  the  for¬ 
ward  simulation  of  Figure  B-3.  It  has  12  outputs,  taken  from 
points  where  the  inputs  were  previously  applied.  This  is  indi¬ 
cated  in  Figure  B-4  with  M  name,  where  M  stands  for  "miss  due  to" 
and  name  indicates  he  source  causing  the  miss. 

The  influence  of  target  maneuvers  and  initial  conditions  on 
miss  distance  is  obtained  directly  after  scaling  at  the  appropri¬ 
ate  outputs  of  the  adjoint  simulation.  The  influence  of  noise 
disturbances  on  miss  is  obtained  by  first  squaring  the  output, 
then  scaling  and  taking  the  integral  of  the  squared  output.  In 
order  to  obtain  the  rms  miss  distance  the  square  root  is  then 
taken.  Thus  for  the  rms  miss  distance  due  to  glint  noise 
obtained  at  A1  in  Figure  B-4,  the  output  is  first  squared  and 
then  integrated  and  scaled,  after  which  the  square  root  is  taken. 

Thus  a  single  run  of  the  adjoint  simulation  provides  the 
effects  on  missile  miss  distance  of  each  of  the  input  quantities, 
such  as  target  maneuvers,  initial  conditions  and  statistical  dis¬ 
turbances,  so  that  one  may  see  which  factors  are  the  most  impor¬ 
tant.  This  represents  a  tremendous  saving  in  ti-me  over  Monte 
Carlo  techniques  (repeated  simulation  trials  plus  ensemble 
averaging) . 
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Figure  B-4:  Adjoint  Missile  Simulation 


